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ABSTRACT 

An  optimal  method  of  filtering  radar  data  is  analyzed  and  experimented. 

First  the  theoretical  bases  of  Kalman's  approach  to  filtering 
problems  is  presented.  Then  the  suitability  of  its  conclusions  to 
the  filtering  of  radar  data  is  studied.  The  form  for  the  filter  needed 
in  such  a  case  Is  derived  and,  finally,  a  simulation  is  made  attempting 
to  show  the  real  usefulness  of  Kalman's  filter. 

It  is  concluded  that  the  method  considered  presents  a  very  prom- 
ising set  of  features  which  may  lead  to  advantageous  practical  use. 

The  filtering  process  Is  based  on  the  recurrent  solution  of  matrix 
formulas.  To  accomplish  a  reasonably  fast  computation  of  the  filtered 
values  the  solution  was  programmed  in  the  FORTRAN  language  and  run  in 
a  1604  CDC  computer. 
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I  .   Introduction 

We  will  be  concerned  in  this  work  with  the  filtering  of  radar  data. 

To  our  knowledge  the  solution  of  this  specific  problem,  in  a  major- 
ity of  cases  at  least,  was  never  attempted  with  the  help  of  o  deep  theo- 
retical foundation.  So  far  the  n<x.,fi   Tracker"  seems  to  have  conveyed 
almost  all  the  possibilities  of  practical  processing  of  radar  data. 

With  the  new  approach  to  the  general  problem  of  filtering  announced 
by  Kalman  [l j  a  new  frame  of  practical  processing  stood  up  to  the  engineer 

The  "  oc,|3  Tracker"  has  no  inherent  means  or  need  of  utilizing  the 
known  statistics  of  the  targets  to  be  tracked  in  any  mathematical  or  even 
rational  v/ay. 

The  Kalman  approach  starts  out  by  quantifying  the  known  parameters 
of  the  signal  generating  process.  Broadly  speaking,  the  class  of  prob- 
lems attacked  by  Kalman  consists  of  filtering  and  predicting  the  output 
of  Gauss -Markov  random  processes  when  its  parameters  and  statistics  are 
known.  The  method  described  by  Kalman  being  particularly  suited  to 
machine  computation  was  immediately  explored  in  more  than  one  field 
of  possible  applications. 

One  branch  of  the  work  done  with  the  Kalman  filter  occupied  control 
engineers.  Another  one  gave  results  in  the  tracking  of  artificial 
satel I ites . 

Ranch  £2]  and  Magi  11  £3  J  followed  Kalman  with  the  analysis  of 
cases  in  which  the  process  parameters  are  random  variables. 

It  will  be  shown  that  Mag  ill's  Technical  Report  is  completely 
suitable  to  provide  a  theoretical  basis  upon  which  a  radar  data  filter 
can  be  built  with  fairly  good  results. 


2,  Statement  of  the  problem. 

The  processing  of  data  provided  by  the  various  types  of  radars 
can  be  readily  implemented  with  the  digital  computer. 

The  processing  device  must  be  able  to  allocate  the  echoes  received 
to  one  of  the  various  targets  being  tracked,  to  decide  on  the  appearance 
of  a  new  target,  to  give  numerical  significance  to  the  distance,  azimuth 
and  other  characteristics  detected  by  the  radar,  to  store  sequentially 
the  information  being  obtained  and  to  provide  all  the  data  stored  to 
any  unit  or  entity  allowed  to  obtain  it. 

This  data  will  be  considered  to  be  present  and  available  for 
f i 1 tering  purposes . 

The  data  being  stored  needs  to  be  handled  in  such  a  way  that  the 
various  tracks  maintained  are  filtered  as  accurately  as  possible.   It 
is  also  considered  that  a  prediction  will  have  to  be  made  and  trans- 
mitted to  some  other  system. 

The  characteristics  of  the  data  provided  by  the  computer  storage 
for  an  eventual  solution  of  the  filtering  problem  are,  obviously,  of 
great  importance. 

Admitting  that  the  allocation  of  echoes  to  targets  was  done 
previously  with  total  confidence,  we  do  not  need  to  consider  but  the 
filtering  of  one  target  track. 

This  way  for  each  target  we  will  have  stored  data  with  the 
following  characteristics, 

a)  Numerical  values  assigned  to  distance,  azimuth  and  altitude. 

b)  Times  at  which  the  position  was  detected. 

c)  If  relevant,  data  on  the  position,  speed  and  course  of  the 
vehicle  where  the  radar  is  mounted. 


Adding  to  this  data  collected  in  the  machine  storage  all  the  char- 
acteristics of  the  target  defining  its  dynamic  model  are  also  necessary 
to  improve  the  tracking  efficiency.  We  may  specify  at  this  point  that 
it  Is  absolutely  necessary  to  have  some  "a  priori"  degree  of  knowledge 
about  the  dynamics  of  possible  targets  when  the  first  point  of  the  track 
is  treated  In  the  filtering  process.  As  will  be  seen,  the  filter  will 
not  start  Its  recurrent  computation  unless  previous  Information  on  the 
characteristics  of  the  possible  targets  Is  available. 

We  need  then: 

d)  A  dynamic  model  capable  of  describing  the  motion  of  the  tracked 
object. 

e)  Some  statistics  of  the  inputs  to  that  model. 

Finally,  the  source  of  the  noise  to  be  filtered  out  has  to  be 
considered.  Noise  will  be  Introduced  by  the  radar  itself,  by  the 
processing  prior  to  the  filter  action,  and  by  any  anomalies  in  the 
propagation  of  the  electromagnetic  waves. 

Let  us  consider  then: 

f)  Statistics  of  the  error  introduced  by  the  equipment. 
With  this  data  available  the  filtering  process  will  aim  at 

reproducing  the  real  track  of  the  target,  eliminating  errors  added 
during  the  acquisition  of  data. 

We  will  see  that  this  objective  Is  not  completely  achieved  but 
that  the  mean  square  error  of  the  positions  will  be  reduced  to  a 
minimum,  with  the  help  of  the  Kalman  filter. 

In  d)  we  referred  to  the  need  of  a  dynamic  model  capable  of  describ- 
ing the  motion  of  the  object  being  tracked.  More  precisely  we  can  say 
that  this  model  Is  going  to  Include  or  use,  in  one  way  or  another,  all 
the  data  gathered  by  us  or  available  in  storage.  We  may  even  say  that 


this  dynamic  model  Is  central  to  the  solution  we  are  seeking.  With 
this  in  view  let  us  proceed  to  formalize  the  model  and  obtain  a  more 
rigorous  concept  of  the  problem  presented  to  us. 

The  signal  or  message  and  the  noise  precesses  are  given  by, 

s(tH)  » <j>t  (t)s(t)  i-  A  (t+i)  U  (t*d 
yd)  -  nYt)s(t) 

w(t«)  ■  <p„  (t)  y(t)  *■  a  (t+i)  uu  (t+i) 

n(t)  ~hr(t)  »(t), 

Where  the  various  symbols  have  the  following  meanings: 

£(f)      Is  the  present  state  vector  of  the  message  process,  i.e., 
a  column  matrix  formed  by  the  target  components  of  distance,  velocity, 
acceleration,  etc.,  on  some  coordinate  axis  (which  we  assume  to  be  the 
X-axis  from  now  on).   In  chapter  three  we  will  Indicate  how  the  other 
components  can  be  treated  by  the  same  filter. 

<fi5(t)     is  the  transition  matrix  of  the  message  process,  interpreted 
as  the  matrix  by  which  we  must  multlpy  the  present  state  to  obtain  the 
next,  as  will  clearly  be  seen  In  the  end  of  this  chapter.  This  includes 
the  knowledge  of  d) . 

Us(t)     Is  a  vector  that  includes  all  the  possible  inputs  to  the 
message  process,  as  wind  gusts  (which  we  will  consider  as  being  able 
to  change  instantly  the  position  of  the  target),  air  speed  and  varia- 
tions of  acceleration  as  caused  by  the  pilot's  action.  As  later  will 
be  seen,  these  inputs  have  to  be  white  gauss  Ian  random  variables  In 
order  that  the  theory  can  be  applied. 

J)s(t)     Is  the  distribution  matrix  of  the  message  process  and 
determines  in  which  way  the  inputs  affect  the  various  state  vector 
components . 


y(t)     Is  the  message  output  which  may  consist  of  one  or  more  of 
the  state  vector  components, 

?(£*■')    is  the  output  vector  of  the  message  process  which  determines 
what  state  vector  components  are  going  to  be  used  in  the  model  output. 

J^v  C%)   's  tne  transition  matrix  of  the  noise  process. 

Uu(£)   is  the  white  noise  process  vector.   Includes  the  knowledge 
of  f). 

ft(fc)   is  the  noise  process  output. 
n(y   Is  the  output  vector  of  the  noise  process. 

The  central  problem  of  this  investigation  can  now  be  clearly  defined, 
The  information  received  by  the  radar  and  sequentially  stored  is  the  sum, 

z(*J~y(t)+n(ti  ,     t -t,*,*, , 

where  /(&)   represents  the  true  distance  and  n(t)   tli.i  irror  added  to 
that.  We  will  seek  a  procedure  of  machine  computation  to  obtain  y(t] 
at  each  instant  £=1,2,3,...  and,  in  a  more  general  fashion,  to  obtain 
some  yftjs y(&)+€y(£)      where  the  error  &/(?/     is  minimized  in  accordance 
with  the  mean-square  error  criterion. 


Notfce  here  that  as  we  are  treating  the  X-axis  components  this 
will  have  to  be  interpreted  as  distance  along  the  X-axis.  This  apparent 
confusion  could  be  avoided  but  we  decided  to  maintain  completely  Hag  ill's 
notation  to  provide  some  unity  between  his  report  and  this  application. 


3.  Consideration  of  the  theory 

Later  on  we  will  study  Kalman  filter's  adequacy  to  radar  data  fitter* 
ing  and  will  concentrate  in  the  areas  of  application.  But,  obviously, 
a  clear  understanding  of  the  theoretical  results  to  be  used  has  to  be 
provided.  This  is  the  present  chapter's  objective. 

The  elements  required  for  a  solution  may  be  concentrated  in  the 
four  equations  presented  in  the  last  chapter  which  describe  the  message 
and  noise  processes. 

A  simplification  of  this  can  be  made  and  consists  of  joining  to- 
gether the  message  and  noise  processes,  reducing  these  to  a  single 
process  which  has  as  output  the  sum  ~z(t)=  yft~)+n(£)    This  is  accom- 
plished by  augmenting  the  state  vector. 

According  to  this  we  now  define: 

Xflt/wf -_/  -  state  vector  of  the  observable  process. 


transition  matrix  of  the  observable  process. 


\  w(t)J 

FU  {  o  Wo) 

**(-W 

hr>@-[r(t)\h(t)j  -  output  vector  of  the  observable  process. 

Defining  also, 
&$vYtfj£Q&J    -  covariance  matrix  of  the  inputs. 


The  mode!  for  the  process  which  generates  the  observed  signal 
may  now  be  reduced  to  the  set  of  equations: 

xfivjm  <t>(t)x(t)  +  *(*+,)  u(tn) 

where      !/(&+?)  ~U(t)  ,   O^T<f 

Even  9f  this  group  of  definitions  casts  some  shadow  on  the  physical 

meaning  of  the  parameters  we  will  not  try  to  clarify  it  right  now.  The 
next  chapter  will  provide  full  understanding  by  means  of  an  example. 
We  may  now  represent  the  whole  process  in  a  block-diagram  form, 
Fig.  1„  where  the  parallel  lines  represent  flow  of  vectors  instead  of 

a  scalar  value  as  in  the  line  at  right  for  2-(&). 


u(ud 


§D(t+/) 


H0^ 


X(t)  M 


*C*J    J 


*»w 


Fig.  1  Hodel  of  the  observable  process 
*As  Q(tJ    Ss  a  diagonal  matrix,  P     (£J    is  equivalent  to  apply 


the  exponent  — L     to  all  non-zero  elements  of  this  matrix, 


Let  us  proceed  now  to  took  at  the  solution  of  the  filtering  problem. 

The  objective  is  to  find  a  recurrence  formula  to  obtain  for  each 
sampling  point  the  best  estimate  of  the  state  vector,  given  both  the 
previous  one  and  observed  state  vectors  at  that  sampling  point. 

The  best  estimate  of  the  state  vector  we  need  is  by  definition, 

k(t/t-i)*£{x(*J/Zt-,}, 

Where  Z t,t     represents  all  the  available  data  up  to  and  including  time  t-l. 

2(t/t-0±  z(i)-z(t/t-,), 

where  z(6)   is  the  observed  signal  at  time  t  and  z(£/t-/)   is  the  best 
estimate  of  the  observed  signal  at  time  t    given  the  previous  best 
estimate.   It  is  possible  to  derive  [3J  a  recurrence  formula, 

x  (Ht)=x  (t/t-0  «  kfct)  2(t/t-i)  ,  t*/,2,s,..., ,     0) 

where  K(t,t)  represents  a  variable  gain  vector. 

We  now  observe  that  the  filter  produces  the  best  estimate  at  time  t 
given  the  data  collected  up  to  and  including  that  time  by  means  of  the 
following  operations, 

a)  Evaluates  a  best  estimate  at  time  t  given  data  collected  up  to 
and  including  t-i  .  This  Is  simply  done  by  multiplying  the  transition 
matrix  <f>(t-l)   by  the  best  estimate  Xft-//t-/J    . 

b)  Corrects  k(£/?-/J    by  the  addition  to  It  of  a  term  containing 
the  information  received  at  t   .  This  term  is  the  product  of  a  weighting 
coefficient  -  the  gain  K(^ttJ   -  and  the  error  In  estimating  the  observed 
signal  given  data  collected  up  to  and  including  t-/   . 

From  what  has  been  said  we  understand  that  the  filter  Is  essen- 
tially formed  with  a  model  of  the  message  process  and  the  added  feature 
of  being  able  to  correct  the  updating  of  Its  output  with  the  information 
received. 
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The  recurrence  formula  (1)  Is  the  solution  to  the  filtering  problem. 
Of  course  the  means  of  updating  the  gain  matrix  is  still  not  uncovered. 

The  derivation  of  a  recurrence  formula  for  this  matrix  starts  out 
considering  that  the  best  estimate  can  be  obtained  by  minimizing  the 
trace  of  the  following  matrix, 

P(t/t)  i  £  {x  (tit J  S  T(t/tjj 

where      X  (tit)  i  X  (t)  -  X  (t/tj  , 

X  (t)   being  the  true  state  vector  of  the  message  and  noise  gener- 
ated. 

The  matrix  P(t/t)   is  the  covariance  matrix  of  the  estimation  error 
and,  obviously,  gives  an  indication  of  the  error  value  at  any  t~  . 

Further  defining, 

r  *  (Ht-i)  *  m  r(tJ  />(&£-/) *,(&)*  m*/2(6/Hj      (2) 

It  is  possible  to  obtain 

k(t,  t)  -  P(t/t-/)  m  (tj  (3) 

A  recurrence  formula  for  Plfc/t-t)  must  still  be  obtained.  Again 
Magi  1 1  derives  it  based  on  the  assumptions  previously  made  and  reaches 
the  result, 

P(tH/t)=t>(t)[T-k(t,  t)rr>Vtj]F>(t/t-t)$r(thl>(9l>TM        ft 

The  recurrent  procedure  is  now  completely  obtainable,  provided  the 
initial  covariance  matrix  P(t9  /  previous  knowledge)  Is  given.  This  is 
done  by  sequentially  evaluating  ^ftjtJ^CtJt), Pft+/J&)        »  etc-« 

With  the  equations  we  have  introduced,  a  filter  block-diagram  is 
possible. 


*(t/t-$ 


mT(t)    t 


j>(t-t)     <fc 


F5g.  2  Model  of  the  optimal  filter 
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4.  Applications 

In  this  chapter  we  intend  to  analyze  more  closely  the  assumptions 
made  before  getting  the  results  of  the  last  chapter  and  test  them  in 
the  problem  of  filtering  radar  data.  We  will  then  decide  on  the 
applicability  to  our  problem  of  the  theoretical  results  presented  before. 

As  It  was  said,  all  Magi  IPs  derivations  are  based  on  the  fact  that 
the  processes  involved  are  either  Gauss-Markov  or  can  be  obtained  from 
these  by  application  of  a  linear  matrix  multiplier.  Can  the  state 
vector  of  any  flying  target,  whatever  dynamic  properties  it  may  have, 
be  considered  a  Gauss-Markov  process?   In  fact,  given  any  instantaneous 
state  vector,  It  Is  always  possible  to  consider  some  other  previous 
time  when  the  definition  applies,  i.e., 

f>[x(**)/x(*J,  *(*-*),-]-  t>[x(t+0/x(*Kx(*-0,^(M)], 

where  K    Is  finite. 

In  a  physical  sense  we  may  say  that  if  we  want,  at  any  Instant  t  , 
to  obtain  the  state  vector  at  t+i    ,  we  may  discard  the  knowledge  we  have 
of  the  target's  state  vector  before  the  instant  t-A:  .  For  example,  if 
we  consider  an  aircraft  at  t   ,  clearly  It  Is  seen  that  the  knowledge 
of  all  its  positions  since  the  flight  started  is  normally  irrelevant  to 
determine  the  next  one  but  we  are  almost  certainly  interested  in  the 
last  two  or  three.  Of  course  the  number  /C    will  vary  with  the  type  of 
target  considered  but  in  any  case  Is  a  finite  number. 

In  Chapter  two  we  said  in  d)  that  we  would  need  a  dynamic  model 
for  the  tracked  target.  This  model  and  the  information  of  e)  (some 
statistics  of  the  inputs)  and  f)  (statistics  of  the  radar  error)  were 
all  considered  in  the  formulas  then  presented. 
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let°s  analyze  first  the  message  process , 

s  (in)  .  <pt  (t)  s  (t) *  4  (tn)  us  ((h) 
y(t)  .  rr(t)s(t) 

If  we  consider  in  the  state  vector  only  61  stance ,  speed  and 
accel Herat  Son  and  this  way  eliminate  the  possibility  of  a  third  deri- 
vative in  the  motion  we  have, 


*„(*)  *Jt)  tJO 


W«)  »J*»)  iM 


Simplifying  the  notation  and  multiplying  matrices  we  obtain, 

The  first  three  terms  of  the  second  member  of  the  equation  repre- 
sent the  formula  to  obtain  the  distance  at  some  instant  £°H    ,  the  dis- 
tance, velocity  and  acceleration  at  t   being  known.  The  components^, 
<pti    and  yn    could  be  easily  gotten  through  the  Taylor  series  expansion 

and  as  we  know  they  are 

where  T~  is  obviously  the  sampling  period. 

The  other  components  of  the  fi   matrix  are  easily  obtained  in  the 
same  fashion. 

We  can  now  give  full  significance  to  QCty/t&C&l)      If   w«  recall 
the  definition  of  the  Ug(pHj     matrix  given  before.   It  is  simply  the 
way  of  correcting  the  next  value  of  the  target's  state  vector  by  adding 
the  possible  random  inputs.  For  example,  we  may  obtain  the  target's 
next  velocity  value  if  we  know  the  present  one,  the  acceleration  and 
the  sampling  period  but  notice  that  we  still  have  to  consider  the  Intro- 
ductlon,  during  the  sampling  period,  of  some  wind  velocity  and  then  we 
have  to  complete  the  value  of  S^Ct+JJ- St(ty+  7 S3CO   by  the  addition 
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of  (JLOth)    . 


The  random  inputs  ^(^*y   an<*  ^aC^O      account  for  instantaneous 
changes  of  position  and  for  variations  of  acceleration  introduced  by 
differences  in  the  engine's  thrust. 

All  put  together  in  the  matrix  equation  we  have, 


s,  (M) 

I        T 

ir* 

s.&) 

/      o     0 

"s,  &*>) 

£*  (tn) 

O       / 

r 

s*W 

o    /   a 

Uu  (tn) 

e/M 

0      o 

/ 

S^Ct) 

o    o    / 

i 

ua(t*0 

One  feature  of  the  random  inputs  when  considered  together  is  that 
one  of  them  may  substitute  the  entire  vector.  Consider,  for  example, 
{/ssC^O  •     Through  this  acceleration  input  we  can  account  for  any 
variation  in  position  or  velocity  which  instantaneously  may  happen. 

Looking  now  at  the  second  equation  of  the  message  process  and 
considering  \/(tj    as  a  vector  containing  only  distance, 


{Kt)] 


o 
o 


In  the  noise  process  we  have  a  different  case.  We  consider  the 
noise  at  any  sampling  time  to  be  Independent  from  any  other  noise  value. 

This  makes  all  the  elements  in  <f>„  (t)   null. 

Let  us  make  it  clear  also  that  the  state  vector  W(t+tJ   has  only 
one  component  which  is,  obviously,  the  error  introduced  in  the  measure- 
ment of  distance,  the  only  one  made.  The  input  (/»(£+*)     's  then  the 
actual  error  V/(t*-/J    . 

We  have  then, 

[»(*»)]-  [o] [w(t)~j  +[,]  [u„  (*,)],  (f) 


and, 


bet)]*  iq[*it}} 


& 
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A  point  we  still  did  not  look  at  is  the  gaussian  quality  of  the 
random   inputs.  Cam  we  admit  that  the  acceleration  imparted  to  an 
aircraft  when  considered  in  discrete  amounts  constitutes  a  random 
variable  with  gaussian  distribution  and  zero  mean?   If  we  think  of 
all  the  aircraft  and  all  its  flights  we,  naturally,  would  be  ready 
to  accept  that.  And  actually,  this  is  what  happens  once  the  filtering 
must  be  appropriate  for  any  aircraft  that  can  appear  in  the  radar  scope, 
This  has  to  be  kept  in  mind  when  values  are  going  to  be  given  to  the 
inputs,  We  will  return  to  this  problem  in  the  next  chapter. 

Having  formed  the  message  and  noise  processes  we  may  now  proceed 
to  construct  the  model  for  the  observable  process,  as  explained  in 
Chapter  2. 

So, 
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and  with  simplified  notation, 


kCtf 


a* 


<4,"*T 


%,<&* 


**J4*      q2 14» 


uaus^ 


u*%,      *»**     <f^     USw 

As  the  random  variables   VU>V**>>  &*s>Us# 
see  that, 

and, 

since  Ust^Us^      »  etc.  are  uncorrelated. 


have  zero  mean,  we 


ffl* 


I  0  0  o 

o  /  o  o 

O  6  /  O 

O  o  0  / 


fauj'*         o  o  o 


o 


f%,  OOO 

o  vv5I  o  o 
o  o  GQ  o 
o        o       o     <n^w 

where  (T  represents   the  standard  deviation. 

With  the  matrix  just  obtained  for  P(0    we  are  able  to  use  the 
recurrence  formulas  presented   in  the  proceeding  chapter. 

But  we  still    need  the   initial    covariance  matrix  without  which  we 
cannot  start  the  computation. 

It  can  be  observed  very  easily  that  at   t\   ,    i.e.,  when  the  first 
measurement  of  the  observable  process   is  made,   the  correlations  between 
distance  and  speed  errors,  distance  and  acceleration  errors,  etc.,  are 
alt    null.     We  have  then, 


P(t*//»*e*''***s  know/eeffe)  = 


*&Kt  O          0  o 

O  4*rZz         o  O 

O  O  **r**  0 

0  0            0  wXj 
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Of  course  the  diagonal  values  must  be  based  on  the  information 

about  the  target  which  we  have  even  before  we  begin  tracking  it.  To 

avoid  a  purely  conjectural  evaluation  of  the  variance  we  may  simply 

discard  the  three  first  received  positions  from  the  radar ,  as  far  as 

filtering  is  concerned,,  and  make  use  of  them  to  estimate  the  velocity 

and  acceleration  of  the  target.  At  t,  +2.    we  will  have  then  the  position, 

the  velocity  and   the  acceleration  of  the  target  with  errors  the  variance 

of  which  we  can  determine. 

Let  us  see  what  are  those  variances.  Supposing  that  we  know  the 

variance  of  the  error  in  the  distance  measurement? 

By  definition  X,  (t+Z)=  Kt(t*+Z)-Xt  (t*+t) 

Then  &r  K,  (t*+L)  ~   l4*r  [K  fa+2J  -  X,  (t.  +2)] 

If  we  take  X*U**2)    as  the  observed  value,  then  we  have, 

&r  X,  (t>+Z)   ->  tfir  U;„ 

because  the  error  in  the  distance  measurement  is  the  error  introduced 

by  the  radar. 

We  want  now  /d/*  Xz(t*+2j   .  Using  the  last  two  points  we  have, 

using  the  same  assumptions  as  above, 

X2  (t*+Z)  -   X,fa*2)-X,{t.*-/)  _    X,(t.  +  2)-X,ft*i-/)    _ 

7"  7" 

-.  x.Ct+z)  -  k.C&+z)  -  rx,(tb+/)-x,C&+*)l  - 

_xtr&+-z)-x,r*,*-/) 

r 

Finally      VdnX2(t*+2)mJ-     &*lx*(h+z)*Xi(k+OJ   "» 

7-2 

Let  us  obtain  now 


The  three  points  we  have  permit  us  to  write, 
j?3  (tM)  =  X5(tQ+2)~  X3  (to+2)  - 

m  x?fa+z)  -  k*CAh)-  *?  Q&lLzJL  fen)  - 

T  T  4  ' 

K(&+z)-x,fa*<)    x,(M-X,fo)  %A#)-x/M  Xfcv)-xM 

T  _  T T  ~  T  _ 

T  T 

-  X/^J-^^^/M^J         X.fa+-2)-ZXtfti>H)+Xt  fa) 


\4tr  £,  (toi-z)m  -L   fen£(h+Z)+4**rX,('&i-/J+iArX,  (h)J  = 


Now 
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5.  Computet  J  on  of  the  filtering  process 

In  the  beginning  of  this  work  we  showed  our  Interest  In  the 

eventual  machine  computation  of  the  recurrent  process,   lit  Is  now 
time  to  look  at  the  procedure  introduced  in  the  last  chapters  and  try 
to  evaluate  how  effective  the  machine  computation  can  be  made, 
in  accordance  with  what  we  mentioned  in  chapter  three  the 
formulas  below  constitute  the  recurrent  process, 

k(t/t)  .  Pfl-/twJ™c& 

The  program  for  this  computation  can  be  done  very  easily.  We 
have  in  the  next  page  a  very  simple  flow-chart  that  illustrates 
this  array  of  computations,, 
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Fig.  3  Flow-chart  for  recurrence  formulas 
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The  problem  of  how  to  start  the  recurrent  process  has  still  to 
be  considered.  We  said  before  that  we  could  obtain  the  value  of  the 
initial  covariance  matrix  by  simply  taking  the  three  first  points  of 
the  raw  track,  evaluating  with  them  the  position,  velocity  and  accel- 
eration at  the  third  point  and,  at  the  sane  time,  calculate  the  initial 
covariance  matrix  of  the  state  vector  thus  obtained. 

The  formulas  then  mentioned  are  repeated  here  for  convenience, 

T 

j    ■'.    ,      vet-  +2) -yfh*)     WkHj-M.) 


7" 


\/+rX,  fa**)*  V**  £/s~ 
tffr  Kg  (tc+2)  m  &•    v&Lh, 


Considering  the  computations  above  we  obtain  the  flow  diagram 
on  the  next  page 
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This  viay  we  use  the  initial  covariance  matrix  at  t^+2     ( JTs  y) 
to  obtain  the  next  one  at  Ta+3     and  use  the  observed  positions  at  U*-2 
to  obtain  a  best  estimate  at  "&i-3  Our  filter  will  then  begin  filtering 

at  and  including  <£>3  . 

To  avoid  too  long  statements  which  would  add  only  to  confusion  in 
the  program  we  divided  the  computation  of  the  covariance  matrix  in 
parts  as  it  can  easily  be  seen. 

Another  point  in  need  of  consideration  is  the  way  the  program  is 
Implemented  without  the  help  of  sub-routines,  even  if  some  matrix 
operations  are  repeated  several  times.  We  thought  that  some  computing 
time  would  be  saved  by  eliminating  the  call  of  sub-routines.  This  way 
there  is  no  great  compilation  time  spent  and  probably  the  compilation 
done  comes  out  more  efficient. 
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6.   Example  of  a  simplified  simulation 

In  this  chapter  we  are  finally  ready  to  implement  the  techniques 
which  have  been  studied  until  now.   Ideally  that  should  be  done  in  a 
real  case.  A  seemingly  feasible  way  would  be  to  obtain  radar  data  from 
maneuvering  airerafts,  project  the  points  in  the  X,Y    and  "2  axis  and 
try  to  filter  the  three  tracks  thus  obtained.  This  way  of  implementing 
our  methods  is,  however,  not  too  flexible. 

A  better  way  is  to  create  a  track,  add  some  noise  and  try  to 
filter  this  out  afterwards.  This  approach  will  provide  us  with  the 
possibility  of  varying  the  track  at  our  will  and  thus  permit  the 
evaluation  of  the  filter's  suitability  in  various  conditions.  How 
much  the  results  can  be  relied  upon  depends,  obviously,  on  the  simu- 
lation's fidelity.  This  we  will  try  to  assess  in  the  next  chapter. 

As  we  have  said  before,  to  study  our  filtering  problem  it  is 
irrelevant  to  consider  more  than  the  projection  of  the  target's  track 
on  one  of  the  coordinate  axis.  The  procedure  for  the  other  axis  is 
exactly  the  same. 

Limiting  ourselves  to  the  projection  on  the  X~axSs  we  can  have 
three  different  situations  -  constant  velocity,  constant  acceleration 
and  sinusoidal ly  varying  distance.  These  correspond  in  three  dimensions 
to  constant  velocity,  constant  acceleration  and  circular  track.  Any 
other  track  of  irregular  form  can  be  reduced  to  a  group  of  circles  with 
accuracy  as  high  as  we  want. 

It  seems  then  only  reasonable  that  we  study  tracks  where  these 
three  possibilities  happen  in  different  orders  and  with  different 
values.   From  here  we  may  eventually  conclude  that  the  filter  provides 
acceptable  results  Sn  any  conceivable  simulated  situation. 
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Furthermore  St  seems  logically  acceptable  to  conclude  the  following: 

a)  The  filter  will  yield  good  results  with  any  constant  velocity 
given  W   that  is  the  case  for  the  lowest  and  highest  expected  velocities 
simulated. 

b)  The  filter  will  work  in  good  conditions  for  any  constant 
acceleration  given  that  ot  works  for  the  highest  possible  and  fulfills 
the  requisite  a) . 

c)  The  filter  will  definitely  accept  any  simulation  and  yield  good 
results  Sf,  furthermore,  it  gives  acceptable  results  with  the  tightest 
circular  maneuver  possible  -  in  the  sense  that  the  acceleration  on  the 
X-axis  reaches  the  highest  value  possible  with  the  speeds  available. 

At  this  point  we  will  have  to  decide  what  types  of  targets  are 
we  going  to  be  handling.  Let  us  confine  ourselves  to  airplanes. 

To  simulate  a  track  we  will  suppose  an  aircraft  going --through  all 
the  situations  presented  above. 

Furthermore  it  seems  desirable  to  have  all  the  combinations  of 
those  pieces  of  track  that  seem  to  introduce  different  degrees  of 
difficulty. 

The  tightest  maneuver  we  want  to  explore  is  the  one  which  gives 
the  greatest  maximum  acceleration  of  the  motion's  projection  on  the 
X-axis,   it  is  easy  to  see  by  Fig.  2-29  in  one  of  the  references  [J* J 
that  the  highest  speed  and  greatest  bank  angle  will  give  the  largest 
value  for  acceleration.  Considering  that  the  maximum  possible  man- 
euvering speed  for  the  aircraft  is  1000  knots  and  the  banking  angle 
70  degrees,  we  have  that  the  projection  on  the  X-axis  of  a  circle 
described  in  the  X~/  plane  Is, 

where  ft  is  the  radius  of  turn  in  miles,,  ^oc*"  the  rate  of  turn  in  degrees 

2k 


per  se<i 

Extrap©lat5n9  Jim  the  tables  for  1000  ml   we  have, 

R-3Z.000    ft-    g^gjg     ft  =6.4  m/ 

fcot9  T.d  x  -^    =.  .  OS237    r*c/ /*ec 

Them  the  maximum  acceleration  is, 

*!/      =  £-(&>t)*~  s^xfoszsy)*- »?;/  sec*  - 

Let   us  assume  that   in  a  straight   line  the  maximum  acceleration 

possible   is   1.805  g°s.     This  gives,, 

/t305  X3Z*  #/sec2=  &**    /»//cec*  m  &S3 x/oJ/w  /sec* 

6K/03  ' 

To  have  the  largest  range  of  speeds  we  will  consider  that  our 
aircraft  can  fly  from  10  to  1000  mi/hrD  which  is  unreal  but  will  give 
greater  scope  to  our  results. 

Let  us  see  now  what  types  of  tracks  will  be  worthwhile  to  exper- 
iment with.  Below  are  the  different  combinations  which  seem  to  present 
different  degrees  of  difficulty  to  our  filter. 

a)  Lowest  speed,,  acceleration  to  highest  speed,,  highest  speed, 
circle  at  highest  speed. 

b)  Lowest  speed,,  acceleration  to  highest  speed,  circle  at  highest 
speed,  highest  speed. 

c)  Inverse  sequence  of  a). 

d)  Inverse  sequence  of  b). 

Another  factor  that  may  be  important  is  the  time  during  which  the 
aircraft  will  stay  at  the  lowest  and  highest  speeds. 

Two  problems  are  still  not  solved.  One  5s  the  addition  of  noise 
to  the  track  we  are  going  to  create.  The  other  is  the  sampling  interval 
which  corresponds  to  the  radar  rotation  period. 
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The  noise  will  be  considered  to  have  a  Gaussian  distribution  with 
a  variance  of  four  miles  constant  for  any  distance.  Of  course  this  is 
a  rather  strong  simplification  and  in  reality  the  variance  will  depend 
very  much  on  the  distance.  We  do  not  think  that  it  is  too  important 
to  deal  with  this  problem  right  now.   In  the  next  chapter  this  will 
be  discussed  again. 

The  sampling  period  will  vary  between  some  limits  In  a  real 
situation.  But,  again,  that  fact  does  not  seem  to  be  important  in 
this  simulation. 

The  programming  of  the  radar  data  simulation  is  very  clear  and 
can  be  analyzed  in  the  Appendix. 

First,  a  time  array  is  generated,  TZMSQc)*     Then,  for  each  piece 
of  track  we  proceed  as  follows: 

For  each  sampling  point  the  noise  is  generated  with  the  RNDEV 
sub-routine  which  provides  us  with  random  numbers  obeying  to  a 
gaussian  distribution  with  variance  of  unity.  DEV  is  the  output 
number,  which  we  multiply  by  VCOE  to  get  a  variance  of  valve  COE. 
Then  we  evaluate  the  signal  output  vector  using  formulas  of  elementary 
physics.  Finally,  the  distance  is  added  to  the  noise  giving  the  input 
to  our  filter. 

A  last  question  is  what  values  to  give  to  the  Q  matrix.  First 
of  all  we  are  going  to  assume  that  all  the  random  inputs  are  null 
except  the  acceleration.  Then  we  wl 1 1  try,  by  the  use  of  several  values 
forQjS  ,  to  obtain  one  which  fits  best  the  whole  filtering  for  all 
the  cases  we  are  going  to  experiment  with. 

The  values  that  finally  were  reached  are  printed  out  In  the  Appendix, 
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To  help  in  the  evaluation  of  the  filtering  quality  a  smoothness 
Index  was  defined  as, 

smooth   ,  *rC**r) ->/c*r$QE)  ^^J 
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7.  Conclusions  and  acknowledgements 

The  results  obtained  indicate  that  the  Kalman  filter  can  be  used 
In  the  filtering  of  radar  data  with  some  success. 

The  cases  taken  as  examples  did  indicate  that  in  a  wide  range  of 
real  situations  the  filter  will  provide  useful  filtering,  as  far  as 
the  simplifications  made  are  not  the  reason  why  the  filter  worked. 

The  constant  sampling  period  which  was  used  could  have  been  varied 
along  the  track  without  any  difficulty.  To  prove  this  we  only  note  that 
for  each  computation  made  after  the  arrival  of  a  new  radar  position  the 
value  used  for  7~8s  not  related  in  any  way  to  the  last  value  of  7" . 

The  radar  error  may  have  a  bias  depending  on  the  distance.  This 
bias  has  to  be  experimentally  determined  and  its  value  subtracted  from 
the  range  measured.   Df  not  it  will  eventually  introduce  a  bias  in  the 
ff 1 tered  values. 

The  Gaussian  quality  of  the  radar  error  is  a  point  of  doubt.  Only 
in  the  presence  ©f  a  definite  radar  system  can  that  be  asserted.  How- 
ever in  a  general  way  we  observe  that  all  components  of  the  radar  noise 
seem  to  have  a  Gaussian  quality. 

The  value  we  should  give  to  the  Q    matrix  is  critical.   Odeally, 
for  each  type  of  target  (he!  3 copters „  airplanes,  missiles,  ships,,  etc.) 
we  should  have  a  Q     matrix.  The  values  for  Q33    can  be  chosen  experi= 
mentally.  The  filter  will  work  with  less  error  when  the  variance  of 
random  inputs  is  small,  as  St  would  be  the  case  of  a  missile  when  com= 
pared  with  an  aircraft,  especially  if  we  consider  a  ballistic  missile. 

The  components  of  the  motion  on  the  three  coordinate  axis  can 
obviously  be  evaluated  in  succession  and  with  a  normal  speed  computer 
the  computation  will  be  completed  before  the  next  set  of  data  arrives. 
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APPENDIX 
OPTIMAL   FILTER   PROGRAM 


PROGRAM    FILTER    3 

DIMENSION    Y(300),G1(300),G2(300),G3(300) , G4 ( 330 ) , YERROR ( 300 ) , 
2X1 (300) , X2(300) ,X3(300) ,XN( 300) ,T IME ( 300 ) , I  TITLE ( 12) • PI  1 (300) , 
3P12( 300), PI  3(300) , P 1 h ( 300 ) , P2 1 ( 300) ,  P22( 300 ) , P23 < 300 ) , P2U ( 300 ) , 
4P3"i  (300),P32(300),P3  3(  3C0  ) ,  P3k  (  300)  ,  PU  1  (  300  )  ,Pu2(300)  ,P43(300). 
5P4M300)  ,XNT(300),X1T(3CO),X2T(300)  , X  1  5 ( 300  )  , X2E( 300)  ,X1ES(3CC)  t 
6AVX1E(500) ,VAXlES(300), VAX! E    ( 3  00 ) , AVXNT ( 300 ) , VAXNT (300 ) , 
7XNTS (300), VAXNT S( 300), SMOOTH ( 300) ,YERRS(300) , A VYERR ( 300 ) , 
8VAYERRS ( 300 )  , VAYERR( 300  ) , ABSXNT ( 300 ) , ABSX IE ( 300 ) 

C  THE    NEXT    CARD    READS    COE    WHICH    IS    THE    SQUARE    ROOT    OF    THE 

C  MEASUREMENT    NOISE    VARIANCE. 

READ    709.C0E 
709    FORMAT(FIO.O) 

READ    2,  P12(3)  ,P13(3),P1M3), 

2  P2K3),  P23(3),P2U{3), 

3  P31 (3),P32(3),  P3U(3), 
k                   Phi (3),P42(3) ,PU3(3) 

2  FORMATM2F5.0) 

READ    3,Q1  l,Q22,C33,Qi+U 

3  FORMAT { h F20.0) 

C  THE    NEXT    CARD    READS    THE    SAMPLING    PERIOD. 

READ    5,T 
5    FORMAT     (F10.0) 

C  THE    NEXT    CARD    READS    THE     INDEXES    WHICH    LIMIT    THE    VARIOUS 

C  PIECES    OF    TRACK. 

READ    300, IMAXl ,  IMIN2, IMAX2,  IM  IN  3, IMAX3, IMIN4, IMAX4 
300    FORMAT     (7110) 

C  THE    NEXT    CARD    READS     IN    SUCCESSION    THE    INITIAL    DISTANCE 

C  AND    SPEED    AND    THE    ACCELLERAT IONS    IN    THE    FOUR    PIECES    OF    TRACK. 

READ    30  3,X1T0,X2T0,X3T0  1,X3T02,X3T0  3,X3T04 
303    F0RMAT(6F10.0) 

C  THE    NEXT    CARD    READS    THE    RADIUS    OF    ROTATION    AND    THE    RATE 

C  OF    TURf4    FOR    THE    CIRCULAR    PIECE    CF    TRACK. 

READ    705, R, ROT 
705    FORMAT     (2F10.0) 

PRINT    999  Cont< 


AH 


999  FORMAT  ( 1  Hi ) 

PRINT  708 
708  F0RMAT(/////10H       CO E  // ) 

PRINT  707, COE 
707  F0RMATIF10.2) 

PRINT  800 
800  FORMAT  (/////7H       T//I 

PRINT  36, T 
36  FORMAT  (F10.2) 

PRINT  301 

301  FORMAT  (/////71H       IMAXl      IMIN2      IMAX2      IMIN3      IKAX 
13      IMINU      IMAXU//) 

PRINT  302, IMAX1.IMIN2, IMAX2 ,  IMIN3.I MAX3, IMIN4, IMAX4 

302  FORMAT  (7110) 
PRINT  30U 

304  FORMAT  (/////72H       X 1T0         X2T0        X3TC1        X3TC2 
2    X3T03        X3T0U   //) 

PRINT  305, XI  TO. X2 TO, X3TC 1 , X3T02 , X3T03, X3T0U 

305  FORMAT  (6F12.8) 
PRINT  306 

306  FORMAT  {/////28H         R  ROT    //) 
PRINT  307, R, ROT 

307  FORMAT  (2F1U.8) 
PRINT  706 

706  FORMAT  (/////56H         Qll  Q22  Q33 

2CUU    //) 
PRINT  31, Qll, Q22,Q33,Q4U 
31  FORMAT  (4F14.8) 
DO  88  K= 1,250 
TIME  (K)=K 
88  TIME(K)=T»TIME(K) 


Cont. 
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C  THE  NEXT  GROUP  OF  C/>RDS  EVALUATES  THE  OUTPUT  OF  THE 

C      OBSERVABLE  PROCESS. X1T  AND  X2T  ARE  THE  REAL  VALUES  OF 
C      DISTANCE  AND  SPEED, XNT  IS  THE  MEASUREMENT  NOISE  AND  Y  IS 
C      THE  OUTPUT  OF  THE  OBSERVABLE  PROCESS. 

NUNIF=1220703125 

DO  105  I=1,IMAX1 

CALL  RNDEV  (NUNIF,CEV) 

XNT ( I)=COE*DEV 

X1T(I)=X1T0       +R*SINF(ROT«(TIME( I  )-TIME(  1     ))) 

X2T(I)=X2T0       «COSF(ROT*(TIME( IJ-TIMECl      ))) 
105  Y( I )=X1T( I)+XNT( I  ) 

DO  102  I=IMIN2, IMAX2 

CALL  RKDEV  (NUNIF.DEV) 

XNT{ I )=COE*DEV 

X1T{ I )  =  X1T{ IMAXl  )  + (TIME (  I  )-TIME(  IMAXl  )  )*X2T(  IMAXl )+ 
2        0.5»(TIME(I)-TIME(  IMAXl)  )*«2*X3T02 

X2T{ I)=X2T( IMAXl )  +  ( T I  ME ( I )-T IME ( I MAX  1 ) >«X3T02 

102  Y( I)=X1T( I)+XNT( I) 
DO  103  I  =  IMIN3t  IMAX3 
CALL  RNDEV(NUNIF,DEV) 
XNT( I )=COE*DEV 

XI T( I)=X1T( IMAX2)+(TIME(I )-TIME( I MAX 2) )*X2T( IMAX2)+ 
2        0.5*(TIME(I)-TIME( IMAX2)  )**2*X3T03 
X2T( I)=X2T( IMAX2)+(TIME( I)- TIME ( IMAX2 ) ) *X3T03 

103  Y( I)=X1T( I )+XNT( I) 
DO  104  I=IMINU, IMAX4 
CALL  RNDEV  (NUNIF,DEV) 
XNTt  I)=COE*DEV 

XI T( I )=X1T(IMAX3)+(TIME(I )-TI ME ( IMAX3 ) )*X2T( IMAX3J+ 
2        0.5«<TIME( I)-TIME( IMAX3)  )**2*X3T04 
X2K  I)=X2T(  IMAX3)+(TIME(I)-TIME(IMAX3)  )»X3T04 

104  Y(  I)=X1T(  I)+XNT(I) 
XI (3)=Y(3) 
X2(3)=(Y(3)-Y(2) )/T 

X3{3)=(((Y(3)-Y{2))/T)-((Y(2)-Y<in/T))/T 
XN(3)=C0 

P11(3)=Q44 

P22(3)=2.0/(T«*2)*G.4U 
P33{3)=6.0/(T»*4)*Q44 
P44(3)=QU4 
X1ES(2)=0.0 
VAX1ES(2)=0.0 
XNTS(2)=0.0 
VAXNTS(2)=0.0 
YERRS(2)=0.0 
VAYERRS(2)=0.0 
600    CONTINUE 
1  =  3 
J=I-1 
GO    TO    10 
21    J=I-1  Cont. 
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All 
Al^ 
A13 
AIM 
A21 
A22 
A23 
A24 
A31 
A32 
A33 
A34 
A41 
A42 
A43 
A44 
Bl  1 
B12 
B13 
B14 
B21 
B22 
B23 
B24 
B31 
B32 
B33 
B34 
BUI 
B42 
B43 
B44 
Cll 
C12 
C13 
C14 
C21 
C22 
C23 
C24 
C31 
C32 
C33 
C34 
C41 
C42 
C43 
C44 
Pll 
P12 
P13 
P14 
P21 
P22 
P23 
P24 
P31 
P32 
P33 
P34 
P41 
P42 
P43 
P44 


:{ 


THE  NEX 
1.0-G1 
•  G 
iO.O 
=-Gl  (J) 
-G2 
1  .0 
=  0.0 

=  -G3(  J) 
=  0.0 

=  1  .c 


T  GROUP 
(J)  ) 


OF  CARDS  EVALUATES  THE  COVARIANCE  MATRIX 


(J) 


U) 


=-G3(J) 
=-G4( J) 
=  0.0 
=  0.0 
=1.0-G4( 
=  P11  (  J)  + 
=P12(J)+ 
=  P  1 3  {  J  ) 
=  0.0 

=P2i (J)+ 
=P21(J)+ 
=P23(J) 
=  0.0 

=P31 { J)+ 
=P32(J)+ 
=P33( J) 
=  0.0 

=  P41 (J)+ 

=P42(J)+ 
=P43(J) 
=  0.0 

=All*Bll 
=A1 1»B12 
=A11*B13 
=  0.0 

=A21*B1 1 
=A21*B12 
=A21*B13 
=  0.0 

=A31*B11 
=A31*B12 
=A31*B13 
=  0.C 

=  A41 *B1  1 
A41*B12 
A41»B13 
O.G 

I)=C1 1+ 
I)=C12+ 
I )=C13+ 
I)=0.0 
I )  =  C2  1  + 
I)=C22+ 
I )=C23+ 
I)=0.0 
I)=C31 
I)=C32 
I )=C33+ 
I)=0.0 
I)=0.0 
I)=0.0 
I )=0.0 
I)=Q44 


J) 

P12( J)*T+P13{J)*0.5»T»»2 

P13( J)*T 


P22( J)«T+P23(J)*0.5*T«*2 
P23( J)*T 


P32( J)«T+P33(J)*0.5*T*»2 
P33(J)*T 


P42(J)*T+P43(J )*0.5*T*«2 
P43{ J)«T 


+  A  14*8*41 
+A14*B42 
+A14*B43 

*B41 

,*B42 

A24*B43 

+A33*B31+A34«B41 
+A33*B32+A34*B42 
+A33*B33+A34*B43 

+A44*B41 
+A44*B42 
+A44*B43 

T*C21+0.5*T**2*C3l4-Qll 

T*C22+0.5*T**2*C32 

T*C23+0.5*T**2*C33 

T*C31 

T*C32+Q22 

T«C33 


Q33 


Cont 
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C  THE    NEXT    GROUP    OF    C4RDS    EVALUATES    THE    GAINS    AND    THE 

C  FILTERED    STATE    VECTOR. 

10    SIGMAS=P1  1( I  )  +  PlU(  I)  +P4K  I  )+PU4(I) 

Gld)  =  (Plld)  +  P14(I))/S  IGMAS 

G2( I )=(P21 (I)+P2M I ) )/S IGMAS 

G3d  )  =  {P31  (  I  )+P34(  I)  )/S  IGMAS 

GU{  I)  =  (P4H  I  )  +  PUU(  I)  )/S  IGMAS 

IF    (  1-3)    200,200,201 

200  YERRORd  )=Y(  I )  -  C  XI  (  I  )+XN(  I)  ) 
XI  (  I  )=X1  (I)+YERRCR(I )*G1(  I) 
X2(  I)  =  X2( I )+YERROR(  I)*G2d  > 
X3(I  )=X3{ I)+YERROR( I )*G3( I) 
XNt  I  )=XNd)+YERROR(  I  )*GM  I) 
GO    TO    202 

201  Y ERROR {  I  ) =Y ( I ) - (X 1 ( J ) +T *X2{ J ) +0 .5«T««2* X3( J ) ) 

XI  ("I  )=X1  (  J)+T»X2(J  )+0.5*T«*2«X3(J)+YERR0R(  I  )  *G  1  ( I  ) 
X2( I  )=X2( J)+T*X3( J  )+YERROR(  I)«G2( I) 
X3(I )=X3( J)+YERROR(I ) *G3( I) 
XN(I )=XN( J)+YERROR( I )*Gk( I) 

202  X1E( I ) =X1 (I)-X1T( I) 
X2E(  I  )=X2(I)-X2T{ I ) 
T0TI=I-2 

XlESd  )  =  X1ES(  J)+X1E(  I  ) 
AVX1EU  )=X1ES(  I  )/TGTI 
VAXlESd  )=VAX1ES(  J  )  +  XTE(I)**2 
VAXlEf I)=VAX1ES ( I )/TOTI 
XNTS( I )=XNTS(J)+XNT( I ) 
AVXNTd  )=XNTS(  I  J/TOTI 
VAXNTSf I)=VAXNTS( J)+XNT (I)**2 
VAXNTt I )=VAXNTS(I)/TOTI 
YERRS( I )=YERRS( J)+YERROR(I) 
AVYERRd  )=YERRS(  I  ) /TOT I 
VAYERRS( I )=VAYERRS(J)+YERROR( I ) *»2 
VAYERR( I)=VAYERRS( I)/TOTI 

SMOOTH ( I )  =  { (VAXNT(  I)-VAX1E(  I) )/VAXNT( I) )»100.0 
ABSXNTd  )=ABSF(XNT(  I  )  1  +  6.0 
ABSXlEt  I)=-ABSF(X1E(  I) )+6.0 
1  =  1  +  1 

IF  (  I-IMAX4-1)  21,22,22  c-. 

22  CONTINUE 
PRINT  999 


Cent . 
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THE  NEXT  GROUP  OF  CARDS  PRINTS  OUT  THE  RESULTS, 
PRINT  756 

756  FORMAT  (//////61H 
\       2  033///) 
1    PRINT  757, Q33 

757  FORMAT(F68.15) 
PRINT  751 

751  FORMAT  (///119H      TIME      XNT        X1T         XI        X1E 

2    X2T         X2        X2E         X3       AVXNT      AVX1E     SMCOT 
3H  //) 

PRINT  750, ( (TIME( I ) , XNT ( I  ) , X1T(  I),X1(I),X1E(I) , X2T( I ) , X2 ( I )  , 
2  X2E(  I),X3( I), AVXNT (I  ) , AVX 1 E ( I ) , SMOOTH! I ) ), I=4,IMAX4) 

750  FORMAT  ( F9. 1 , 1 1  El  0.3 ) 

PRINT  999 

PRINT  752 

752  FORMAT  (//////69H      TIME      Pll        P22        P33        Gl 
2      G2         G3    //) 

PRINT  7  53,((TIME(I),P11  m,P22<  I )  ,P33(  I)  ,G1  ( I  )  ,G2(  I )  ,G3(  I  )  )  , 
2  I=4.IMAX4,2) 

753  FORMAT  (F9.1  ,6E10.3) 
DO  1 1  L=l  ,12 

11  ITITLE(L)=8H 

ITITLEd  )=8HCARL0S  P 
ITITLE(2)=8H.  SIMOES 
ITITLE(4)=8H    MS   T 
ITITLE(5)=8HHESIS 
ITITLE(7)=8H       OP 
ITITLE(8)=8HTIMAL  FI 
ITITLE(9)=8HLTERING 
ITITLE(10)=8H0F  RADAR 
ITITLEt 1 1 )=8H  DATA 
LABEL-MH  Pll 

CALL  DRAW  ( 1 49, T I  ME, P 1 1, 1,0, LABEL,  I  TITLE, 0, 0,0,0,0,0,5,6,1,  LAST  I 
LABEL=4H  P22 

CALL  DRAW  (  1  49 , TIME, P22, 2,0 ,L ABEL ) 
LABEL=UH  P33 

CALL  DRAW  ( 1 49, T IME, P33, 3,0 , LABEL ) 
LABEL=4HVXNT 

CALL  DRAW  ( 249, TIME, VAXNT , 1 ,0, LABEL, ITITLE , 0,0, 0, 0, 0,0,5,6 , 1 .LAST) 
LABEL=MHVX1E 

CALL  DRAW  (249, TIME, VAX1E, 2,0, LABEL ) 
LABEL=4H 

CALL  DRAW  ( 249 , TIME, AVXNT ,2 ,0, LABEL ) 
LABEL=4H 

CALL  DRAW  ( 249, TI ME, AVX 1 E,3 ,0, L ABEL ) 

CALL  DRAW(249,TIME,ABSXNT, 1,0, LABEL, ITITLE, 0,0 ,0,0, 0, 0,5,6, 1 , LAST) 
LABEL=UHX1E 

CALL  OR AW ( 249, T IME, ABSX IE, 3,0, LABEL ) 
LABEL=4HSM00 

CALL  DRAW (249, TIME, SMOOTH ,0,0, LABEL, ITITLE, 0,0,0,0,0,0,5,6 , 1 ,LAST) 
601  CONTINUE 
END 
END 


A~6 


INPUT  DATA 


COE 
2.00 

T 
2.00 


IMAX1  IMIN2      IMAX2      IMIN3      IMAX3      IMINU      IMAXU 

60  61        148        1 1*9        162        163        250 

X1T0  X2T0        X3T01        X3T02        X3T03        X3T04 

7.00000000  .27800000    .OOCOOOOO    .00000000   -.00983000    .00000000 

R  ROT 

6.U0000C00  .05237000 

Qll  Q22            Q3i            QUH 

.OOOOOOOC  .OOCOOOOO      .0000C900     4.OCG00C0O 
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Case  d) 
OUTPUT  DATA 


TIME  XNT  X1T  Xl  X1E  X2T 


8.0 

.263E+01 

.898E+01 

.106E+02 

. 166E+01 

•264E+00 

10.0 

-.115E-01 

.960E+01 

.971E+0 1 

. 103E+00 

.254E+00 

12.0 

-.158E+00 

.  102E+C2 

.976E+01 

-.4145E+00 

.214  1E+00 

14. 0 

•224E+01 

.108E+02 

. 123E+02 

. 151E+01 

.225E+00 

16.0 

-.2ii7E+01 

.113E+02 

.989E+01 

-. 140E+01 

.207E+00 

18.0 

-.190E+01 

.  118E+02 

.958E+01 

-.218E+01 

.186E+00 

20.0 

-.1 19E+01 

.  122E+02 

. 102E+02 

-.  195E+01 

.163E+00 

22.0 

-.2C9E+01 

.  125E+02 

.  102E  +  02 

-.230E+01 

.139E+00 

2U.0 

-.  1j>6E  +  01 

.128E+02 

.  109E+02 

-. 198E+01 

.1  13E  +  00 

26.0 

-.150E+01 

.131E+02 

.113E+02 

-. 181E+01 

.858E-01 

28.0 

.8ii2E+00 

.133E+02 

.129E+02 

-.406E+00 

.577E-01 

30.0 

-.1hlE+01 

. 134E+02 

.127E+02 

-.687E+00 

.290E-01 

32.0 

.394E+01 

. 134E+02 

.151E+02 

. 174E+01 

-.844E-04 

34.0 

-.335E+00 

. 134E+02 

•148E+02 

. 141E+01 

-.291E-01 

36.0 

.425E+01 

. 133E+02 

.164E+02 

•315E+01 

-.579E-01 

38.0 

-.198E+01 

. 151E+02 

.148E+02 

.175E+01 

-.860E-01 

40.0 

.8^6E+00 

.128E+02 

.147E+02 

. 181E+01 

-.1  13E  +  00 

42.0 

.418E+01 

.125E+02 

.157E+02 

.312E+01 

-.139E+00 

44.0 

-.152E+01 

.122E+02 

.141E+02 

. 195E+01 

-.163E+00 

46.0 

-.190E+01 

. 118E+02 

.126E+02 

.884E+00 

-.186E+00 

48.0 

-.1G6E+01 

. 113E+02 

.  116E+G2 

.332E+00 

-.2C7E+00 

50.0 

-.316E+00 

. 108E+02 

.-109E  +  02 

. 137E+00 

-.225E+00 

52.0 

.215E+01 

. 102E+02 

.110E+02 

.783E+00 

-.24 1 E+00 

54.0 

-.761E+00 

.960E+01 

.993E+01. 

.327E+00 

-.254E+00 

56.0 

-.926E+00 

.897E+01 

.886E+01 

-. 1 15E+00 

-.264E+00 

58.0 

•215E+01 

.833E+01 

.880E+01 

.478E+00 

-.272E+00 

60.0 

-.951E+00 

.767E+01 

.764E+01 

-.291E-01 

-.276E+00 

62.0 

.850E+00 

.700E+01 

.709E+01 

.924E-01 

-.278E+00 

64.0 

.385E-01 

.633E+01 

.625E+01 

-.765E-01 

-.276E+00 

66.0 

.482E+00 

.567E+01 

.557E+01 

-.977E-01 

-.272E+00 

68.0 

-.131E+01 

.502E+01 

.434E+01 

-.679E+00 

-.264E+00 

70.0 

-.222E+01 

.439E+01 

.294E+01 

-.  145E  +  01 

-.254E+00 

72.0 

.413E+00 

.380E+01 

.251E+01 

-. 129E+01 

-.241E+00 

74.0 

-.191E+00 

.323E+01 

.191E+01 

-.  133E+01 

-.225E+00 

76.0 

-.184E-01 

.271E+01 

.141E+01 

-. 131E+01 

-.206E+00 

78.0 

-.624E+00 

•224E+01 

.778E+00 

-.  146E  +  01 

-.186E+00 

80.0 

-.145E+01 

.182E+01 

-.171E-01 

-. 184E+01 

-.163E+00 

82.0 

.700E+00 

. 145E+0  1 

-.214E-01 

-. 148E+01 

-.139E+00 

84.0 

-.145E+00 

.  115E+01 

-.268E+00 

-. 142E+01 

-.1 13E+00 

86.0 

-.  171E+01 

.912E+00 

-.922E+00 

-. 183E+01 

-.857E-01 

88.0 

. 123E+01 

.739E+00 

-.527E+00 

-.127E+01 

-.576E-01 

90.0 

.262E+01 

.634E+00 

.286E+00 

-. 348E+00 

-.288E-01 

92.0 

.215E+00 

.600E+00 

.289E+00 

-.311E+00 

.253E-03 

94.0 

.168E+01 

.636E+00 

.807E+00 

. 171E+00 

.293E-01 

96.0 

.288E+01 

.741E+00 

.168E+01 

.941E+00 

.581E-01 

98.0 

.396E+01 

.915E+00 

.282E+01 

.191E+01 

.862E-01 

A-8 


X2 


X2E 


X3 


AVXNT 


AVX1E 


SMCOTH 


.332E+00 

.703E+CC 

.UUlh+CC 

.662E+00 

.443E+C0 

.U21E+C0 

.155E+00 

. 129E+00 

. 1U6E-01 

.761E-01 

.31UE+00 

.218E+00 

.523E+00 

.379E+C0 

.515E+00 

.239E+00 

.173E+C0 

.246E+C0 

.U08E-C1 

. 128E+00 

.226E+GC 

•.280E+C0 

.25CE+C0 

.327E+CC 

.395t+CC 

•360E+CC 

.428E+CC 

.430E+C0 

•.455E+CC 

•.461E  +  CC 

■.514E  +  00 

.575E+C0 

-.539E+C0 

•.520E  +  C0 

•.491E  +  C0 

•.475E  +  00 

■.1474E  +  C0 

-.4G1E  +  C0 

•.355E  +  00 

-.351E  +  00 

■.251E  +  CG 

■.124E+C0 

-.836E-01 

■.207E-02 

.1G2E+CC 

.217E+00 


.673E-01 

•.957E+00 

•.682E  +  00 

.437E+00 

.649E+00 

•.607E+00 

•.319E  +  00 

-.268E  +  00 

■.984E-01 

•.975E-C2 

.256E+C0 

.189E+00 

.523E+00 

.408E+00 

.573E+C0 

.325E+00 

.286E+00 

.385E+00 

.2OUE+C0 

.585E-01 

•.  194E-01 

■.552E-01 

•.957E-02 

-.732E-01 

-.  130E+00 

-.882E-01 

-.151E+00 

•.152E+00 

■.  179E+00 

-.190E  +  G0 

■.250E  +  00 

-.321  E+OO 

-.299E+00 

-.295E  +  00 

■.285E+00 

■.289E+00 

-.311  E  +  OO 

■.262E  +  C0 

-.21+2E+00 

■.265E  +  00 

-.194E  +  00 

-.953E-01 

-.839E-01 

-.3UE-01 

.439E-01 

.131 E+GO 


■.  196E  +  CC 

-.360E+C0 

■.  159E  +  00 

.797E-C1 

•.981E-01 

•.717E-01 

•.287E-01 

•.209E-C1 

■.533E-C2 

.  165E-03 

.167E-01 

.8U0E-02 

.258E-01 

.  147E-01 

.204E-01 

.431E-C2 

.653E-03 

.402E-02 

-.572E-02 

-.127E-G1 

-.160E-01 

-.169E-01 

-.142E-01 

-.163E-01 

-.178E-01 

-.H*7E-C1 

-.164E-C1 

-.150E-C1 

■.1U8E-01 

-.138E-01 

-.  1U9E-01 

■.163E-C1 

-.  132E-01 

-.  112E-C1 

■.883E-C2 

-.731E-02 

-.661E-02 

■.267E-02 

-.368E-C3 

■.138E-03 

.U37E-02 

.973E-02 

.107E-01 

.134E-01 

.  169E-01 

.206E-01 


.170E+01 

.113E+01 

.805E+00 

.  109E+01 

.499E+00 

.157E+00 

.114E-01 

.243E+00 

.355E+00 

-.U59E  +  00 

•.350E+00 

■.U40E+00 

.127E+00 

•.H+1E+00 

.134E+00 

.932E-02 

.542E-01 

.271E+00 

.182E+00 

.826E-01 

.305E-01 

.15UE-01 

.104E+00 

.698E-01 

.315E-01 

.1  lOE+OO 

.721E-01 

.989E-01 

.969E-01 

.109E+00 

.650E-01 

-.418E-02 

.809E-02 

.241E-02 

.183E-02 

•.151E-01 

-.528E-01 

-.335E-01 

-.363E-01 

-.771E-01 

■.459E-01 

.162E-01 

.207E-01 

.577E-01 

.119E+00 

.201E+00 


.121E+C1 

.839E+00 

.518E+C0 

.716E+00 

.364E+00 

.607E-03 

.2U3E+00 

.472E+00 

.623E+00 

.730E+00 

.703E+00 

.702E+00 

.527E+00 

■.398E  +  00 

.176E+G0 

.630E-01 

.413E-01 

.203E+00 

.291E+00 

.319E+00 

.320E+00 

.312E+00 

.331E+00 

.331E+00 

.3H4E  +  00 

.320E+C0 

.308E+G0 

.300E+00 

.288E+00 

.275E+C0 

.245E+00 

.19UE+00 

.150E+00 

.108E+G0 

.688E-01 

.275E-01 

.216E-01 

-.589E-01 

•.929E-01 

•.135E+00 

•.162E+00 

•.167E+00 

-.170E+00 

■.162E+00 

-.138E  +  00 

■.948E-01 


.558E+02 
.557E+02 
.532E+02 
.538E+02 
.584E+C2 
.1+38E  +  C2 
.311E+02 
.22  9E+C2 
.146E+C2 
.1C4E+02 
.118E+02 
.162E+C2 
.359E+02 
.321E+C2 
.355E+02 
.3U8E+02 
.3C9E+02 
.335E+02 
.31 1E+02 
.328E+C2 
.335E+02 
.235E+02 
.359E+02 
.362E+C2 
.367E+C2 
.392E+02 
.397E+02 
.U01E+C2 
.UC1E+02 
.U02E+02 
.U07E+02 
.U1UE+02 
.U01E+02 
.386E+02 
.371E+02 
.355E+02 
.338E+02 
.323E+02 
.306E+C2 
.295E+02 
.291E+02 
.327E+02 
.327E+C2 
.341E+02 
.373E+C2 
.U12E+02 


A=9 


TIME      XNT        X1T         XI        X1E        X2T 


100.0 

.274E+01 

. 116E+01 

.348E+01 

.232E+01 

.TT3E+00 

102.0 

.521E+01 

. 146E+01 

.486E+01 

.340E+01 

.139E+00 

104.0 

-.4C3E+01 

.183E+01 

.329E+01 

.  146E+01 

.164E+00 

106.0 

,2!j3E+00 

.225E+01 

.337E+01 

. 1 13E+01 

•186E+00 

108.0 

.205E+00 

.272E+01 

.356E+01 

.834E+00 

.207E+00 

1  10.0 

.171E+01 

.324E+01 

.43OE+0  1 

. 1C6E+01 

.225E+00 

112.0 

-.30  1E+01 

.381E+01 

.362E+01 

-. 187E+00 

.241E+00 

114.0 

-.192E+01 

.440E+01 

.354E+C1 

-. 868E+00 

.254E+00 

116.0 

-.270E+01 

.503E+01 

.339E+01 

-. 164E+01 

.264E+00 

1  18.0 

-.150E+01 

.568E+01 

.381E+01 

-. 187E+01 

.272E+00 

120.0 

-.343E+00 

.634E+01 

.469E+01 

-. 165E+01 

.277E+00 

122.0 

-.228E+01 

.689E+01 

.497E+01 

-. 192E+01 

•277E+00 

124.0 

-.154E+01 

.744E+01 

.556E+01 

-. 189E+01 

.277E+00 

126.0 

.85  6E-01 

.800E+01 

.667E+01 

-. 133E+01 

.277E+00 

128.0 

-.772E+00 

.855E+01 

.745E+01 

-. 1  10E+01 

.277E+00 

130.0 

-.163E+01 

.910E+01 

.795E+01 

-. 1 16E+01 

.277E+00 

132.0 

.982E+00 

.966E+01 

.926E+01 

-.4C0E+00 

.277E+00 

134. 0 

-.398E+01 

.102E+02 

.893E+01 

-. 128E+01 

.277E+00 

136.0 

-.632E+00 

.  108E+02 

.977E+01 

-.992h+00 

.277E+00 

138.0 

-.152E+00 

.113E+02 

.107E+02 

-.605E+00 

.277E+00 

140.0 

-.154E+01 

. 119E+02 

.  112E+02 

-.710E+00 

.277E+00 

142. 0 

-.101E+01 

.124E+02 

.118E+02 

-.648E+00 

.277E+00 

144.0 

-.168E+01 

.  130E  +  02 

.122E+02 

-. 808E+00 

.277E+00 

146.0 

.434E+01 

. 135E+02 

.  144E+02 

.884E+00 

.277E+00 

148.0 

-.188E+01 

.141E+02 

.145E+C2 

.377E+00 

.277E+00 

150.0 

.155E+01 

.146E+02 

. 156E+02 

.969E+00 

.277E+00 

152.0 

.103E+01 

.152E+02 

.165E+C2 

. 128E+01 

.277E+00 

154.0 

.236E+01 

. 157E+02 

. 177E+02 

. 197E+01 

.277E+00 

156.0 

-. 17  1E+01 

.163E+02 

. 175E+02 

. 120E+01 

.277E+00 

158.0 

.3U9E+00 

.168E+02 

. 180E+02 

. 1  15E+01 

.277E+00 

160.0 

-. 166E+01 

. 174E+02 

.179E+02 

.468E+00 

.277E+00 

162.0 

-.2C3E+01 

. 180E+02 

.177E+02 

-.237E+00 

.277E+00 

164.0 

.532E+00 

.185E+02 

. 185E+02 

-.54  1E-01 

.277E+00 

166.0 

.134E+0  1 

.191E+02 

. 194E+02 

. 339E+00 

.277E+00 

168.0 

-.499E-01 

.196E+02 

.198E+02 

.236E+00 

.277E+00 

170.0 

-.144E+C0 

.202E+02 

.203E+02 

. 1  13E  +  00 

.277E+00 

172.0 

-.247E+01 

.207E+02 

.200E+02 

-.7C4E+00 

.277E+00 

174.0 

-.449E+01 

.213E+02 

.193E+02 

-.200E+01 

.277E+00 

176.0 

.143E+01 

.218E+02 

.206E+02 

-. 125E+01 

.277E+00 

178.0 

.257E+01 

.224E+02 

.221E+02 

-.264E+C0 

.277E+00 

180.0 

.250E+01 

.229E+02 

.235E+02 

.551E+00 

.277E+00 

182.0 

-.472E+01 

.235E+02 

.225E+02 

-.974E+00 

.277E+00 

184.0 

.162E+01 

.2i*0E+02 

.237E+02 

-.309E+00 

.277E+00 

186.0 

.109E+01 

.246E+02 

.247E+02 

.893E-01 

.277E+00 

188.0 

-.983E+00 

•251E+02 

.249E+02 

-.208E+00 

.277E+00 

190.0 

-.162E+00 

.257E+02 

.255E+02 

-.206E+C0 

.277E+00 

192.0 

.8U7E+G0 

.262E+02 

.264E+02 

. 105E+00 

.277E+00 

194.0 

-.378E+01 

.268E+02 

.258E+02 

-. 105E+01 

.277E+00 

196.0 

-.163E+00 

.274E+02 

.265E+02 

-.899E+00 

.277E+00 

198.0 

.209E+01 

.279E+02 

.278E+02 

-.776E-01 

.277E+00 

A-10 


X2 


X2E 


X3 


AVXNT 


AVX1E 


SMOOTH 


.275E+00 
.390E+C0 
.221E+C0 
.2168+00 
.217E+00 
.267E+CC 
.181E+0C 
.155E+00 
• 125E+0C 
.148E+C0 
.210E+CC 
.21 1E+00 
.238E+GC 
.309E+00 
.341E+CC 
.3U3E+00 
.418E+CC 
.334E+C0 
.363E+00 
.397E+00 
.381E+CC 
.381E+00 
.359E+C0 
.506E+00 
.439E+00 
.476E+CC 
.483E+00 
.522E+00 
.42HE+00 
.399E+C0 
.319E+00 
•2UUE+00 
.258E+CC 
.290E+C0 
.274E+00 
.259E+00 
.  182E+C0 
.708E-01 
. 156E+00 
.257E+00 
.333E+00 
■.189E+CC 
.258E+00 
.296E+C0 
.268E+CC 
.270E+CC 
.299E+00 
.  193E+00 
.215E+C0 
.296E+00 


.162E+00 

.250E+00 

.572E-01 

.293E-01 

.101 E-01 

.417E-01 

-.601  E-01 

■.989E-01 

-.139E+00 

-.  124E+00 

•.664E-01 

-.659E-01 

-.388  E-01 

.320 E-01 

.648E-01 

.670E-01 

.1U1E+00 

.577E-01 

.861+E-01 

.120E+00 

.ICUE+OO 

.  10UE+00 

.827E-01 

.230E+00 

.163E+00 

.200E+00 

.2C6E+00 

.2U5E+00 

.1U7E+00 

. 123E+00 

.U25E-01 

-.322 E-01 

-.183E-01 

. 131 E-01 

-.246E-02 

-.179  E-01 

-.942  E-01 

-.206E  +  00 

-.120E+00 

-.  198  £-01 

.561 E-01 

-.873E-01 

-.186E-01 

.199 E-01 

••841  E-02 

-.698E-02 

.223E-01 

-.838  E-01 

-.619E-01 

.195E-01 


.213E-01 

.246E-01 

.1U7E-01 

.132E-01 

.120E-01 

.132E-01 

.810E-02 

.621E-02 

.429E-02 

.U9UE-02 

.729E-C2 

.666E-02 

.728E-C2 

.983E-02 

.104E-C1 

.958E-02 

.121 E-C 1 

.720E-02 

.78UE-02 

.866E-02 

.717E-C2 

.651E-02 

.U94E-02 

.  112E-C1 

.711E-02 

.815E-02 

.770E-02 

.877E-02 

.35UE-02 

.210E-02 

•.172E-02 

-.U95E-C2 

-.390E-02 

-.210E-02 

-.261E-02 

■.307E-02 

-.625E-02 

-.1C7E-01 

■.590E-02 

■.8C8E-03 

.271E-02 

-.UOUE-02 

-.557E-03 

.123E-C2 

•.  159E-03 

-.798E-0U 

. 125E-C2 

-.367E-02 

-.23UE-02 

.156E-02 


.254E+00 

.355E+00 

.267E+00 

.267E+00 

.266E+00 

.293E+00 

.232E+00 

.192E+00 

.  1U1E+00 

.  112E+00 

.104E+00 

.637E-01 

.370E-01 

.378E-01 

.2U7E-01 

•.152E-02 

. 139E-01 

-.U77E-01 

-.565E-01 

•.580E-01 

■.797E-01 

-.933E-01 

-.116E  +  00 

•.531E-01 

-.785E-01 

-.562E-01 

■.416E-01 

-.685E-02 

•.293E-01 

-.24UE-01 

•.453E-01 

-.704E-01 

•.629E-01 

-.U56E-01 

-.U56E-01 

-.U68E-01 

•.757E-01 

■.128E+00 

-.109E+00 

■-787E-01 

•.U93E-01 

-.102E  +  00 

-.826E-01 

■.698E-01 

-.797E-01 

-.806E-01 

-.707E-01 

•.llOE+00 

-.llOE  +  00 

.877E-01 


,41+uE-or 

.259E-01 

.5146E-01 

.756E-01 

.902E-01 

.108E+00 

. 103E+00 

.853E-01 

.545E-01 

.208E-01 

.795E-C2 

.U04E-01 

.712E-01 

.919E-01 

.108E+G0 

.125E+00 

.129E+00 

. 1M7E+00 

.160E+00 

.166E+C0 

. 17UE+00 

.181E+00 

.190E+00 

.175E+00 

.167E+00 

•.152E  +  00 

•.132E+00 

•.lOi+E  +  00 

.873E-01 

•.713E-01 

-.643E-C1 

•.665E-01 

.66UE-01 

•.61UE-01 

■.577E-01 

.557E-01 

.63UE-01 

•.862E-01 

-.997E-01 

.102E+00 

•.91+2E-01 

-.lOUE+00 

•.106E  +  00 

•.104E  +  00 

•.105E  +  00 

•.106E+00 

.104E+00 

.11UE+00 

.122E+00 

.122E+00 


TITO  6  E+ 0  2 
.U30E+02 
.U63E+02 
.U58E+02 
.U5UE+C2 
.U57E+C2 
.479E+02 
.48UE+02 
.U89E+C2 
.479E+C2 
.U67E+02 
.i463E  +  02 
.U5UE+02 
.447E+02 
.14U3E  +  02 
.UUME+02 
.UU5E+02 
.U72E+C2 
.U70E+02 
.U68E+C2 
.471E+02 
.U72E+C2 
.U75E+02 
.506E+C2 
.512E+02 
.513E+02 
.509E+02 
.507E+C2 
.507E+02 
.5C3E+C2 
.506E+02 
.513E+C2 
.513E+02 
.516E+02 
.515E+02 
.515E+02 
.523E+02 
.5U0E+02 
.538E+02 
.5U6E+02 
.553E+C2 
.577E+02 
.580E+02 
.581E+02 
.582E+C2 
.582E+02 
.583E+02 
.595E+02 
.593E+C2 
.598E+02 


A-11 


TIME 


XNT 


X1T 


XI 


X1E 


X2T 


200.0 

202.0 

204.0 

206.0 

208.0 

210.0 

212.0 

214.0 

216.0 

218.0 

220.0 

222.0 

224.0 

226.0 

228.0 

230.0 

232.0 

234.0 

236.0 

238.0 

240.0 

242.0 

244.0 

246.0 

248.0 

250.0 

252.0 

254.0 

256.0 

258.0 

260.0 

262.0 

264.0 

266.0 

268.0 

270.0 

272.0 

274.0 

276.0 

278.0 

280.0 

2  82.0 

284.0 

286.0 

288.0 

290.0 

292.0 

294.0 

296.0 

298.0 


-.8015+00 

-.194E+00 

.181E+00 

-.269E+01 

.174E+01 

-. 166E+01 

.197E+01 

.217E+01 

.1 19E+01 

-.252E+00 

-.407E-01 

-.618E+00 

.417E+00 

.218E+01 

.254E+01 

-.357E-01 

.324E+00 

-.771E+00 

-.293E+C1 

.226E-02 

-•137E+00 

-.274E+01 

. 157E+00 

.314E+01 

.295E+01 

-.186E+00 

.427E+00 

.154E+G1 

.113E+01 

.653E+00 

. 1 55E+00 

.120E+00 

.3Q7E+01 

-.613E+00 

-.354E+00 

-.250E+01 

.160E+01 

.146E+01 

.712E+00 

-.185E+01 

-.301E+00 

.254E+01 

-.373E+01 

-.151E+01 

.142E+00 

.7GOE+00 

.326E+01 

-.944E+00 

.262E+00 

-.192E+01 


.285E+02 

.290E+02 

.296E+02 

.301E+02 

.307E+02 

.312E+02 

.318E+02 

.323E+02 

.3295+02 

.334E+02 

.340E+02 

.345E+02 

.351E+02 

.356E+02 

.362E+C2 

.368E+02 

.373E+02 

. 379E+02 

.384E+02 

.390E+02 

.395E+02 

.401E+02 

.406E+02 

.412E+C2 

.417E+02 

.423E+02 

.428E+02 

.434E+C2 

.439E+02 

.445E+02 

.451E+02 

.456E+02 

.462E+02 

.467E+02 

.473E+02 

.478E+02 

.484E+02 

.489E+02 

.4955+02 

.500E+02 

.506E+02 

.511E+02 

.517E+02 

.522E+02 

.528E+02 

.533E+02 

.539E+02 

.545E+02 

.550E+02 

.555E+02 


.282E+C2 

.288E+02 

.295E+02 

.292E+02 

.305E+02 

.306E+02 

.319E+02 

.332E+02 

.340E+02 

.343E+02 

.347E+02 

.349E+02 

.355E+02 

.367E+02 

.378E+02 

.380E+02 

.384E+02 

.384E+02 

.379E+02 

.384E+02 

.390E+02 

.387E+02 

.395E+02 

»412E+02 

.426E+02 

.429E+02 

.435E+02 

.UU4E+0  2 

.451E+02 

.456E+02 

.460E+02 

.463E+02 

.476E+02 

.477E+02 

.479E+02 

.475E+02 

.485E+C2 

.494E+02 

.500E+02 

.498E+02 

.502E+02 

.51 5E+02 

.503E+02 

.510E+02 

.518E+02 

.5272+02 

.544E+02 

.546E+02 

.552E+02 

.552E+C2 


-.269E  +  00 

-.243E  +  00 

-.  106E  +  00 

-.  869E+00 

-.  149E  +  00 

-.581E  +  00 

. 165E+00 

.849E+00 

. 1 11E+01 

.869E+00 

. 715E+00 

.398E+00 

.439E+C0 

. 1C0E+01 

. 157E+01 

. 125E+01 

. 106E+01 

.566E+00 

-.  513E  +  00 

•.  518E+00 

-.5515  +  00 

-.  136E+01 

-.  1  13E+01 

•.977E-02 

.878E+00 

.662E+00 

.661E+00 

.995fc+00 

. 1  14E  +  01 

. 1  10E  +  01 

.910E+00 

.729E+00 

. 147E+01 

.9555+00 

. 594E+00 

-.3655  +  00 

.885E-01 

.433E+00 

. 4975+00 

-.2305  +  00 

■.357E+00 

.4155+00 

-.  846E+00 

•.121E+01 

■.993E  +  00 

.620E+00 

.488E+00 

. 134E+00 

.206E+00 

.3875+00 


.277E+00~ 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.2775+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.277E+00 

.257E+00 


A-12 


X2 


X2E 


X3 


AVXNT 


AVX1E 


SMOOTH 


.278E+00 
.281E+C0 
.294E+00 
.224E+00 
.295E+00 
.255E+00 
.325E+CC 
.383E+00 
.397E+00 
.362E+00 
.339E+CC 
.302E+00 
.301E+00 
.348E+C0 
.390E+CC 
.3U7E+C0 
.320E+00 
.267E+0C 
.167E+CC 
.175E+C0 
.  180E+00 
.1  15E+0C 
.152E+C0 
.267E+00 
.350E+00 
.325E+00 
.321E+C0 
.348E+00 
.354E+00 
.3U3E+C0 
.318E+00 
.296E+00 
.360E+CC 
.303E+00 
.266E+00 
. 178E+00 
.227E+00 
.263E+C0 
.269E+00 
.204E+00 
. 199E+00 
.277E+C0 
.164E+00 
. 143E+00 
.177E+00 
.223E+00 
.332E+GC 
.298E+00 
.30UE+00 
.248S+C0 


.1U5E-02- 

.460E-02 

.174E-01 

•.527E-01 

.18UE-Q1 

■.217E-01 

.487E-01 

.107E+G0 

.120E+00 

.859E-01 

.622E-01 

.255E-01 

.245 E-01 

.710E-01 

.113E+00 

.701E-01 

.437 E-01 

-.908E-02 

-.109E  +  00 

-.102E  +  00 

•.968E-01 

-.161  E+OO 

-.125E  +  00 

•.994E-02 

.734E-01 

.U83E-01 

.446E-01 

.711E-01 

.775 E-01 

.662E-01 

.412E-01 

.19UE-01 

.832E-01 

.268E-01 

-.105E-01 

-.983  E-01 

-.493E-01 

-.140E-01 

-.719E-Q2 

-.726E-01 

•.772E-01 

.792E-03 

-.113E  +  00 

-.134  E+OO 

-.992  E-01 

-.531  E-01 

.555E-01 

.211 E-01 

.279 E-01 

-.872E-02 


1  Jy  /L   U-» 


.688E-C3 

.121E-02 

.208E-02 

.133E-02 

•.606E-03 

.264E-02 

.502E-02 

.517E-02 

.316E-02 

.  180E-02 

•.284E-04 

•.686E-04 

.204E-02 

.378E-02 

.1U8E-02 

.146E-03 

•.226E-02 

-.660E-C2 

■.566E-02 

■.492E-02 

•.739E-02 

-.508E-02 

.593E-03 

.431E-G2 

.279E-C2 

.237E-02 

.335E-02 

.334E-02 

.252E-02 

.117E-02 

.707E-C4 

.295E-02 

.133E-03 

■.157E-02 

■.5U1E-02 

■.269E-02 

■.853E-03 

-.467E-03 

■•339E-02 

-.329E-02 

.542E-03 

•.465E-02 

-.519E-G2 

-.315E-02 

-.771E-03 

.422E-02 

.228E-02 

.238E-02 

■.386E-03 


.960E-01 

.932E-01 

.1 19E+00 

.  101E  +  00 

.116E+00 

.957E-G1 

.741E-01 

.621E-01 

.639E-01 

•.637E-01 

.688E-01 

.644E-01 

.442E-01 

•.21  1  E-01 

.213E-01 

.182E-01 

■.248E-01 

•.498E-01 

.493E-01 

.501E-01 

•.726E-01 

•.707E-01 

•.U42E-01 

•.197E-01 

-.210E-01 

•.174E-01 

•.492E-02 

.U11E-02 

.922E-02 

.10UE-01 

.112E-01 

.348E-01 

.298E-01 

.269E-01 

.792E-02 

.198E-01 

.30UE-01 

.351+  E-01 

.216E-01 

.193E-01 

.375E-01 

.106E-01 

-.251E-03 

.748E-03 

.56UE-02 

.282E-01 

.215E-01 

.232E-01 

.994E-02 


i123E+00 
.125E+00 
.  12UE  +  00 
.132E+00 
.132E+00 
.  136E+00 
.133E+00 
.124E+00 
.112E+00 
.  103E+00 
.956E-01 
.911E-01 
.862E-01 
.764E-01 
•617E-01 
.501E-01 
.404E-01 
•351E-01 
.392E-01 
.433E-01 
.476E-G1 
.586E-01 
.676E-01 
.671E-01 
.593E-01 
.535E-01 
.477E-01 
.394E-01 
.300E-01 
•211 E-01 
.138E-01 
.805E-02 
.332E-02 
.106E-01 
.150E-01 
.121 E-01 
.127E-01 
.158E-01 
.19UE-01 
.176E-01 
.148E-01 
. 177E-01 
.1 15E-01 
.286E-02 
.U15E-02 
.8U6E-02 
.501E-02 
.405E-02 
.261E-02 
.523E-02 


» 596E  +  rs> 
.598E+02 
.598E+C2 
.6C3E+02 
.606E+C2 
.608E+02 
.61  1E  +  C2 
.614E+0  2 
.612E+C2 
.61  1E  +  02 
.61CE+C2 
.610E+02 
.609E+C2 
.61 1E+C2 
.611E+02 
.608E+02 
.605E+C2 
.605E+02 
.612E+02 
.61 1E+02 
.61  1E  +  02 
.613E+02 
.610E+02 
.619E+02 
.62*4E  +  C2 
.623E+02 
.622E+02 
.622E+02 
.62CE+C2 
.618E+C2 
.616E+02 
.615E+C2 
.618E+02 
.617E+C2 
.616E+02 
.621E+C2 
.623E+02 
.62UE+02 
.62UE+02 
.626E+C2 
.626E+02 
.631E+02 
.639E+02 
.638E+C2 
.636E+C2 
.636E+02 
.6U2E+02 
.643E+C2 
.643E+C2 
.645E+02 


A-13 


TIME 


XNT 


X1T 


XI 


X1E 


X2T 


300.0 
302.0 
30*+. 0 
306.0 
308.0 
310.0 
312.0 
314.0 
316.0 
318.0 
320.0 
322.0 
324.0 
326.0 
328.0 
330.0 
332.0 
334.0 
336.0 
338.0 
340.0 
342.0 
344.0 
346.0 
348.0 
350.0 
352.0 
354.0 
356.0 
358.0 
360.0 
362.0 
364.0 
366.0 
368.0 
370.0 
372.0 
374.0 
376.0 
378.0 
380.0 
382.0 
384.0 
386.0 
388.0 
390.0 
392.0 
394.0 
396.0 
398.0 


-.145E+01 

-.706E+00 

•632E+00 

•.808E  +  00 

-.138E+01 

.1G1E+01 

■.14  1E+00 

.540E+00 

.6i8E+00 

.144E+01 

-,2";2E+00 

.985E+00 

■•937E-01 

-.15  3E+01 

.203E+01 

.176E+01 

.3U9E+01 

-.2^5E+01 

-.688E+00 

.126E+01 

-.248E+01 

-.4i)3E+00 

.144E+01 

.260E+00 

.730E+00 

.429E+G1 

.685E+00 

-.330E+01 

-.106E+01 

-.168E+01 

.120E+01 

-.4<6E+00 

.358E+01 

.  175E+C1 

-.175E+01 

-.234E  +  01 

.1C5E+01 

-.114E+01 

.728E+00 

-.223E+01 

-.587E+00 

-.779E+00 

-. 197E-01 

-.135E+01 

.226E+01 

-.3<7E+01 

.974E+C0 

.142E+01 

.275E+01 

.604E+00 


.560E+02 
.565E+02 
.569E+02 
.573E+02 
.576E+02 
.579E+02 
.582E+02 
.584E+02 
.586E+02 
.587E+02 
.588E+02 
.589E+02 
.589E+02 
.589E+02 
.589E+02 
.589E+02 
. 589E+02 
.589E+02 
.589E+02 
.589E+02 
.589E+02 
.569E+02 
.589E+02 
.589E+02 
.589E+02 
.589E+02 
.589E+02 
•589E+02 
.589E+02 
.589E+02 
.589E+C2 
.589E+02 
.589E+02 
.589E+02 
.589E+02 
.590E+C2 
.59CE+02 
.590E+02 
.590E+O2 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 


•553E+02 
.558E+02 
.566E+02 
.569E+02 
.570E+02 
.579E+02 
.582E+02 
.586E+02 
.592E+02 
.598E+02 
.598E+02 
.601E+02 
.599E+02 
.593E+02 
.599E+02 
.602E+02 
.610E+02 
.59  8E+0  2 
.593E+02 
.594E+02 
.58UE+02 
.582E+02 
•586E+02 
.586E+02 
.587E+02 
.600E+02 
.599E+02 
.586E+02 
.582E+02 
.577E+02 
.582E+02 
.581E+02 
.59  3E+02 
.597E+02 
.590E+02 
.583E+02 
.587E+02 
.584E+02 
.587E+02 
.581E+02 
.581E+02 
.580E+02 
.582E+02 
.580E+02 
.589E+02 
.580E+02 
.586E+02 
.592E+02 
.60  1E+0  2 
.601E+02 


.710E+00 
.725E+00 
.303E+00 
. 372E+00 
•.  593E  +  00 
.398E-01 
.775E-01 
. 380E+00 
.661E+00 
.  1  12E  +  C1 
.981E+00 
.  121E+01 
. 105E+01 
.436E+00 
.9712+00 
. 130E+01 
.2C6E+01 
.853E+00 
.350E+00 
.510E+00 
•.5C6E+00 
■.729E  +  00 
-.  330E+00 
-.  344E+00 
•.  1  85E+00 
.  103E  +  01 
.9452+00 
■•  343E+00 
-.731E  +  00 
•.  123E+01 
■.736E+00 
-.804t+00 
.323E+00 
.776E+00 
.857E-01 
-.669E+00 
•.255E  +  00 
•.536E  +  00 
-.222E+00 
-.8532+00 
•.869E  +  00 
■.927E  +  00 
-.7  312+00 
-.962E  +  00 
-.452E-01 
-.928E+00 
-.  400E+00 
.182E+00 
. 107E+01 
. 1  14E+01 


.237E+00 

.218E+00 

.1 98E+00 

.178E+00 

.159E+00 

.139E+00 

.1  19E  +  00 

.996E-01 

.799E-01 

.6C3E-01 

.4C6E-01 

.2C9E-01 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-G2 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

.127E-02 

•127E-02 

.127E-02 


AH4 


X2 


X2E 


X3 


AVXNT 


AVX1E 


SMOOTH 


.218E+00 

.21 5E+00 

.249E+0C 

.232E+C0 

. 199E+00 

.236E+00 

.227E+C0 

.231E+C0 

.230E+00 

.242E+00 

.195E+C0 

.182E+00 

. 132E+0O 

.450E-01 

.710E-01 

.770E-01 

. 123E+00 

•.138E-01 

■.723E-C1 

•.637E-01 

-.  161E+00 

•.175E  +  00 

-.129E  +  G0 

-.123E  +  C0 

•.  101E  +  00 

.165E-C1 

.550E-C2 

-.1 13E+00 

•.  139E+00 

■.  171E+00 

-.109E  +  CC 

•.  102E+C0 

. 123E-C1 

.550E-G1 

-.102E-C1 

■.756E-01 

-.287E-01 

-.489E-01 

•.132E-01 

■.666E-01 

-.589E-C1 

•.557E-C1 

-.295E-01 

-.444E-G1 

.462E-01 

-.353E-C1 

.186E-01 

.721E-01 

. 148E+00 

.141E+CC 


-.192E-01 

■.229E-02 

.511 E-Ol 

.535E-01 

.402E-01 

.969E-01 

.107E+00 

.132E+CO 

.151E+00 

.182E+00 

.154E+00 

.161 E+OO 

.131E+00 

.437 E-Ol 

.698E-01 

.757E-01 

.122 E+OO 

■.  150E-01 

■.736E-01 

-.650  E-Ol 

-.162E+00 

-.176E+00 

-.130E  +  00 

-.124E  +  G0 

-.102E+00 

.152 E-Ol 

.U22E-02 

-.115E  +  00 

-.  140  E+OO 

•.172E+00 

-.1  10E  +  00 

-.104E+00 

.1  10E-01 

.538E-01 

-.115E-G1 

•.768E-01 

-.299 E-Ol 

-.501  E-Ol 

-.145E-01 

-.678E-01 

-.6 02 E-Ol 

-.569E-01 

•.307E-01 

-.457E-01 

.449E-01 

-.365E-01 

.174E-01 

.708 E-Ol 

.146 E+OO 

.140E+00 


.172E-C2 

•.169E-02 

•.362E-05 

•.788E-C3 

•.221E-02 

■.329E-C3 

•.721E-C3 

•.424E-03 

•.439E-03 

. 129E-C3 

■.2C2E-02 

•.242E-02 

■.447E-02 

-.801E-02 

■.610E-02 

•.528E-02 

-.272E-02 

•.866E-02 

•.1C5E-01 

-.919E-02 

-.127E-01 

-.122E-01 

•.904E-02 

-.796E-02 

•.622E-C2 

-.354E-03 

-.821E-03 

-.613E-02 

-.673E-02 

-.756E-02 

-.4C8E-02 

•.31+2E-02 

.209E-02 

.384E-02 

.533E-03 

•.247E-02 

-.126E-C3 

-.103E-02 

.680E-03 

-.180E-02 

-.129E-C2 

•.1C3E-02 

.254E-03 

-.4U5E-03 

.370E-02 

•.326E-03 

.215E-02 

.437E-02 

.739E-02 

.643E-02 


.651E-04 
.U67E-02 
.432E-03 
.578E-02 
.148E-01 
.818E-02 
.904E-02 
.550E-02 
.124E-02 
.791E-02 
.652E-02 
.127E-01 
.  120E-01 
.242E-02 
.149E-01 
.256E-01 
.467E-01 
.316E-01 
.272E-01 
.346E-01 
•196E-01 
.169E-01 
.253E-01 
.266E-01 
.310 E-Ol 
.556E-01 
.593E-01 
.401E-01 
.338 E-Ol 
.241E-01 
.307E-01 
.281 E-Ol 
.467E-01 
.561E-01 
.462E-01 
.332E-01 
.387E-01 
.329E-01 
.366E-01 
.245E-01 
.212E-01 
.170E-01 
.168E-01 
.964E-02 
.214E-01 
.484E-02 
.984E-02 
.171 E-Ol 
.310E-01 
.339E-01 


•.999E-02 

•.148E-01 

•.167E-01 

-.191E-01 

-.228E-01 

•.230E-01 

•.223E-01 

•.197E-01 

-.153E-01 

•.811E-02 

•.185E-02 

.575E-02 

.122E-01 

. 149E-01 

.208E-01 

.286E-01 

.410E-01 

.460E-01 

.1478E-01 

.506E-01 

.472E-01 

.426E-01 

.U05E-01 

.382E-01 

.369E-01 

.426E-01 

.478E-C1 

.456E-01 

.412E-01 

.340E-01 

.297E-C1 

.251E-01 

.267E-01 

.308E-01 

.311 E-Ol 

.273E-01 

.258E-01 

.228E-01 

.214E-01 

.168E-01 

.120E-01 

.708E-02 

.320E-02 

-.186E-02 

-.208E-02 

-.688E-02 

-.891E-02 

-.793E-02 

•.244E-02 

.336E-02 


.645E+02 
.645E+C2 
.645E+02 
•645E+02 
.646E+02 
.646E+C2 
.646E+02 
.646E+02 
.646E+02 
.645E+02 
.643E+02 
.641E+02 
.639E+C2 
.640E+02 
.641E+C2 
.64CE+02 
.640E+C2 
.643E+02 
.643E+Q2 
.643E+02 
.647E+Q2 
.646E+C2 
.647E+02 
.647E+02 
.647E+C2 
.656E+C2 
.655E+02 
.661E+02 
.661E+C2 
.66CE+02 
.660E+02 
.659E+02 
.665E+C2 
.665E+C2 
.667E+02 
.669E+02 
.670E+C2 
.670E+02 
.670E+02 
.671E+02 
.670E+02 
.669E+C2 
.668E+02 
.668E+C2 
.670E+02 
.674E+C2 
.674E+02 
.675E+C2 
.677E+02 
.675E+C2 


A-!5 


TIME 


XNT 


X1T 


XI 


X1E 


X2T 


400.0 
402.0 

uou.o 

406.0 
408.0 
410.0 
412.0 
414.0 
416.0 
418.0 
420.0 
422.0 
424.0 
426.0 
428.0 
430.0 
432.0 
434.0 
436.0 
438.0 
440.0 
442.0 
444.0 
446.0 
448.0 
450.0 
452.0 
45U.0 
456.0 
458.0 
460.0 
462.0 
464.0 
466.0 
468.0 
470.0 
472.0 
474.0 
476.0 
478.0 
480.0 
482.0 
484.0 
486.0 
488.0 
490.0 
492.0 
494.0 
496.0 
498.0 
500.0 


.1955+01 

.265E+01 

.991E+00 

.255E+01 

.266E+01 

.126E+01 

.413E+01 

.91 1E+00 

.175E+01 

.284E+00 

.139E+01 

.26UE+01 

.1 10E+C1 

.166E+01 

■.8'j2E  +  00 

•.1S9E  +  01 

•108E+01 

.103E+0  1 

.231E+00 

•.328E  +  01 

.580E-01 

•.156E  +  00 

-.142E+00 

■.203E  +  00 

.l^ilE+01 

•319E+01 

-.159E-01 

■.334E+01 

.133E+01 

-.186E+00 

.59  1E+00 

■.551E+00 

-.739E+00 

.150E+00 

.143E+01 

.186E+00 

■.214E+01 

.178E+01 

-.490E+01 

.190E+01 

-.426E+00 

.2^8E+01 

-.179E+01 

-.383E+00 

-.321E+Q0 

.234E+01 

. 116E+01 

■.636E+00 

•189E+01 

-.139E+01 

.106E+01 


.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.5  90E+0  2 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
.590E+02 
•591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
•591E+02 
•591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 
.591E+02 


.594E+02 
.602E+02 
.597E+02 
.588E+02 
.596E+02 
.599E+02 
.585E+02 
.588E+02 
.594E+02 
.59  4E+0  2 
.597E+02 
.604E+02 
.604E+02 
.606E+02 
.600E+02 
.593E+02 
.595E+02 
.596E+02 
.595E+02 
.583E+02 
.584E+02 
.584E+0  2 
.58UE+02 
.58UE+02 
.589E+02 
.579E+02 
.581E+02 
.573E+02 
.580E+02 
•582E+02 
.586E+02 
.586E+02 
.585E+02 
.588E+02 
.593E+02 
.594E+02 
.588E+02 
.595E+02 
.580E+02 
.588E+02 
.588E+02 
.597E+02 
.59  1E+02 
.59  0E+0  2 
.590E+02 
.597E+02 
.600E+02 
.597E+02 
.602E+02 
.596E+02 
.598E+02 


.403E+00 

. 1  17E+01 

.690E+00 

-.2014E+00 

. 632E+00 

.896E+00 

•.537E+00 

■.200E+00 

.3U1E+00 

.3U6E+00 

.683E+00 

.  13UE+01 

. 140E+01 

.  161E  +  01 

. 100E+01 

.255E+00 

.440E+00 

.581E+00 

.U55E+00 

-.721E+00 

-.674E  +  00 

-.677E+00 

-.665E  +  00 

-.656E+00 

-.  1U2E  +  00 

-.  1  10E  +  01 

-.926E+00 

-.  177E+01 

-.  103E  +  01 

-.859E  +  00 

-.463E+00 

••473E+00 

-.535E+00 

-.318E  +  00 

.262E+00 

.360E+00 

-.275E  +  00 

.381E+00 

-.  1  11E+01 

-.286E+00 

-.303E+00 

.567E+00 

-.  171E-01 

-.855E-01 

-.  128E  +  00 

.6U2E+00 

.918E+00 

.585E+00 

. 106E+01 

.446E+00 

.666E+00 


.127E-02 
.127E-02 
.127E-02 
•127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
•127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 
.127E-02 


A°16 


X2  X2E         X3  AVXNT  AVX 1 E  SMCOTH 

.6G6E-01  .593E-G1       .220E-02  .239E-01  -"8E-02  .677E  +  02 

.124E+00  .122E+00       .486E-C2  .37TE-0  .     12E-0  .678E+C2 

.665E-01  .653E-01       .183E-02  .320E-01  . U6E-01  -678E+C2 

-.230E-G1  -.2U2E-01    -.239E-02  .191E-01  .135E-0]  .681E+02 

fP8E-0  I515E-01       .126E-02  .322E-01  .166E-01  .68UE+02 

.697E-0  .684E-0    .  91E-02  .382E-01  .209E-01  -683E+C2 

-.691E-C  -I704E-01    -.455E-02  .178E-01  .182E-01  .690E+C2 

-.54    E-Cl  -I  351+E-OI    -.255E-02  .221E-01  .171E-01  .690E+02 

168E-C1  .155E-0      -.167E-04  .305E-01  .187E-01  .692E+02 

I    43E-0  I    30 E-0      -.    27E-C3  I318E-01  .2C3E-01  -691E+C2 

.U23E-0  iJlOE-O         .     15E-02  .383E-01  .235E-0]  -692E+02 

964E-01  I951E-01       .350E-02  .508E-01  .298E-01  .692E+02 

914E-01  .902E-01       I296E-02  .558E-01  .363E-01  .690E+02 

I993E-0  980E-0            305E-02  .634E-01  .438E-01  .688E+C2 

I321E-0  I303E-01  -.276E-03  .591E-01  .483E-0  .687E+C2 

-  419E-C1  -4^2 E-0 1  -.360E-02  .514E-01  .493E-01  .688E+02 

-  239E-01  -.252E-0  -.2S6E-02  I561E-01  .5UE-01  .688E+02 
-IllOE-O  -  123E-0   -I165E-C2  .607E-01  .536E-01  -688E+C2 

-  232E-C  -I2U5E-01  -.206E-02  .615E-01  .554E-01  .688E+02 
-'.129E  +  00  -:i30E  +  00  -667E-02  .U6QE-01  .519E-01  .69  £402 
-.113E+00  -.115E+00  -.535E-02  .461E-01  .485E-01  .69TE+C2 
-I  03E+00  -IlOUE+OO  -.U38E-C2  .U53E-01  .£52E-0  -690E+02 

-  906F-01  -.919E-01  -.344E-02  .444E-01  .420E-01  .69CE+C2 

-  794E-G1  -  807E-01  -I262E-02  .433E-01  .388E-01  .689E+02 

-  230E-0  -.243E-0  -17UE-03  .U95E-01  .380E-01  -690E+C2 
-.105E+00  -.107E+00  -.358E-02  .349E-01  -329E-01  -692E+C2 

-  7A^F-N  -  777F-01  -.194E-02  .347E-01  .286E-01  .691E+02 
:;ISll+00  -IlJuE+OO  -  476E-02  .flll-O  .206E-01  .692EJ02 
-I582E-01  -.595E-01  -.506E-03  .255E-01  .159E-01  .69TE+C2 

-  525E-0  -I338E-0  .704E-03  .246E-01  .121E-01  .69CE+C2 
107E-C  I947E-02  I260E-02  .271E-01  .100E-01  .69CE+C2 
128E-01  IlUE-Ol   I2U6E-02  .245E-01  .789E-02  .690E+02 

.966E-02  .838E-02   .209E-G2  .2J2E-01  -553E-02  -690E+C2 

^TAP-m  ^O^F-Ol   .290E-02  .217E-01  .413E-02  .6S9E+C2 

•.839E-0  I826E-0    .5C0E-02  .278E-01  .524E-02  .690E+02 

869E-01  .857E-Q1   .469E-02  .284E-01  -676E-02  .690E+02 

.224E-01  .211 E-0 1   .134E-02  .192E-01  .556E-02  .692E+C2 

ISOUE-O  792E-0    .385E-02  .267E-01  .716E-02  -693E+C2 

-.622E-01  -I635E-0   -.296E-C2  .579E-02  .2U1E-02  .700E+02 

I188E-01  .175E-01   .975E-03  . 1 38 E-0 1  .120E-02  .701E+02 

I  59E-0  I  U6E-0    .754E-03  .120E-01  -.77UE-04  .701E+02 

I93UE-01  I921E-01   .U20E-02  .223E-01  .230E-02  .7C3E+C2 

.31UE-C  30  E-0    .101E-02  .U7E-01  .222E-C2  -70UE+02 

A    6E-01  .203E-01   .U73E-03  .131 E-0 1  . 1 85E-02  -704E+02 

I  1 U9E-01  1 136E-01   .126E-C3  .117E-01  .131E-02  .70UE  +  C2 

.825E-0  I812E-0    .318E-C2  .213E-01  .395E-02  -706E+C2 

.985E-C  .972E-0    .362E-02  .260E-01  .770E-02  -705E+02 

.573E-0  .560E-01   .1U2E-02  .233E-01  .101E-01  -705E+02 

I929E-0  I916E-01   .290E-02  .308E-01  .143E-0  'l^E*02 

I258E-01  I2U6E-01  -.395E-03  .251E-01  .  161E-01  .7C5E+02 

.406E-01  I393E-01   .310E-03  .292E-01  .187E-01  .705E+02 


Case  d) 

P  MATRICES 


^ 

1 

§ 

CM 

cs> 
(5 

CS) 
09 

I 

Pll 

Si 

&\ 

PM 

MPl                      WH                      flf»                      BM                     ,«* 

X-SCflLE  -  1.00E+62  UNITS/INCH. 
Y-SCALE  -  L00E+81  UNITS/INCH. 

CARLOS  P»  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-18 


Averages  and   variances 

of 

measyrement  noSse  amid  filter  error 


J/WT 


\M1E 


Avxvr 


WXIE 


•  "M"carv5W* 


09i 


:, -SCALE  -  1.08E+82  UNITS'INCH. 
Y -SCALE  -  L88E+80  UNITS/INCH. 

CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


AH  9 


Case 


Absolute  values 
of 

mofse  measurement  amd  filter  erroo 


XHT 


X-SCflLE  -  L08E+02  UNITS'INCH. 
Y-SCflLE  -  2.00E-f8a  UNITS'JNCK 

CARLOS  Po  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A=20 


Case  d) 
Per forma nee 


ioorw 


X -SCALE  -  1.80E+8Q  UNITS'INCH. 
Y-SCflLE  -  Q.88E+81  UNITS'INCH. 

CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING' OF  RADAR  DATA 


A 


Case  b) 


rn 


CSJ 


B 


\ 

ll 

I 

I  V 

pi: 

DM 

f\V 

oee 


mi 


08Q 


883 


084 


X-SCflLE  -  1.88E+82  UNITS'IUCH. 
Y-SCPLE  -  J.88E+81  UNITS' INCH. 

CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-22 


\M\T 


Alt 


X-3CfiLE  -  L08E+82  IHITf/INCH. 
Y-SCflLE  -  1.88E+88  UNITS'INCH. 


CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-23 


Case  b) 


X -SCALE  -  1.00E+Q2  UNITVJtICH. 
T-?CflLE  -  2.8GE4Q0  LMTVIMCH. 


CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-24 


Case 


<r-t 


to 


=H 


<a 


;    1 

9 

\ 

1 

• 

. 

r 

» 

- 

riOOTH 


oee 


081 


oqq 


Q80 


OQ1! 


X-SCflLE  -  L08E+02  UNIT5/INCH. 
Y-*CflLE  -  i.eeE+01  UNITS'IrtCH. 


CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-25 


Case  a) 


g 


CM 

cs 
o 


cs 

Q 


C5> 

5 

cs 


; 

- 

-   j 

I 

Pli 

PHI 

f\V 

088 


681 


08Q 


eeo 


084 


X-SCflLE  -   108E-t82  LUITVINCH. 
Y-SCfiLE  -  L08E481  UNITVINCH. 


CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-26 


g 


J/MHT 


;KI£ 


,AVXlE 


X-SCALE  -  1.O8E+02  UNITS'INCH 
Y -SCALE  -  J.68E-f88  UNITS'INCH. 

CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-27 


,as<g  a 


sKNT 


X "SCALE  -  1.80E+82  UNITS'INCH. 

y-scale  -  2.ee£H-ee  units/'Inch. 

CARLOS  P.  SIMOES  MS     THESIS 

'  OPTIMAL  FILTERING  OF  RADAR  DATA 


Case  a' 


CSJ 


1 

1 

i 

"***"   ^  ■  i,r^ 

i 

1 

MGOfW 


aaa 


88  J 


882 


880 


884 


X -SCALE  -  1.08E+QQ  UNITS' JMCH. 
Y-SCflLE  -  Q.88E+8J  UNITS'INCH 


CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-29 


la%e  c) 


n 


S 


^ 

: 

- 

ft 

A 

pii 

paa 

oae 


081 


002 


080 


884 


X-SCflLE  -  1.88E+8Q  UNITS'INCH. 
Y-SCflLE  -  1.88E+81  UNITS'INCH. 

CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RflDflR  DATA 


Case  c) 


.UXHT 


j;xie 


X -SCALE  -  L88E-f8Q  HUTS'INCH 
Y-SCflLE  -  l.Q8E+d8  UNITS/INCH 


CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-31 


>ase 


X -SCALE  -  i.eeE-t-82  UNITJVJNCH. 

'Y-scflLE  -  ^.eoE-tee  units/inch. 

CARLOS  P.  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A°32 


Case  c) 


OS 


• 

! 

^JJ 

fi 

1 



5MQOTH 


888 


881 


G8Q 


080 


884 


X-3CALE  -   i.08E-fQ2  UNITVINCH 
Y-SCflLE  -  2.88E+81  UNITS'INCU 

CARLOS  Po  SIMOES  MS     THESIS 

OPTIMAL  FILTERING  OF  RADAR  DATA 


A-33 


thesS4945 

Optimal  filtering  of  radar  data. 
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